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I.  INTRODUCTION 


The  objective  of  this  research  program  is  to  develop  a  cumulative 
damage  model  which  can  be  used  to  relate  the  strength,  stiffness,  and 
service  life  properties  of  epoxy-based  structural  composite  laminates.  The 
desired  cumulative  damage  model  would  take  into  account  the  effect  of  load- 
history,  characteristics  of  the  applied  stress  (tension,  compression,  etc.),  lam¬ 
ination  stacking  sequence,  the  effect  of  natural  and  man-made  flaws  and  the 
inherent  mechanical  properties  of  the  basic  composite  material  system.  In 
essence,  a  unified  methodology  for  laminate  failure  analysis  is  to  be  de¬ 
veloped  in  this  research  program. 

The  basic  concept  adapted  in  the  model  development  is  to  regard  the 
accumulation  of  damages  in  composite  laminates  as  a  load  and/or  time  depen¬ 
dent  "process."  The  process  involves  a  certain  type  of  sub-laminate  cracks 
which  initiate  at  locations  where  some  dominating  material  flaws  exist. 

Under  sustained  loading  these  sub-laminate  cracks  would  grow  and  some  would 
coalesce  to  form  a  large  laminate  crack.  At  this  stage,  failure  of  the  lami¬ 
nate  can  be  expected. 

The  term  "sub-laminate  cracks"  pertains  to  small  matrix  cracks  which 
generally  occur  internal  to  the  laminate.  The  origin  of  sub-laminate  cracks 
is  thought  to  stem  from  sub-critical  growth  and  coalescence  of  inherent 
material  microflaws  which  exist  throughout  the  laminate.  These  microflaws 
would  act  as  local  stress  risers  when  external  loading  is  applied  or  a  tem¬ 
perature  change  is  incurred.  At  a  certain  location  where  the  stress  condition 
reaches  a  critical  value,  some  dominant  flaws  begin  to  coalesce  and  form  a 
large  plane-crack.  The  size  of  such  a  crack  is  orders  of  magnitude  larger 
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than  tliat  of  the  microflaws.  It  is  usually  detectable  at  the  macroscopic 

scale,  e.g.,  by  x-radiography ,  sonic  scan  and  even  visual  inspection  fl].  > 

Because  the  matrix  material  and  the  matrix-fiber  interface  are  weaker 
compared  to  the  strength  of  the  fibers,  a  sub-laminate  crack  is  almost 
always  propagated  through  the  matrix  phase  of  the  material.  Hence  it  is 
commonly  labeled  as  a  matrix  crack. 

Since  laminates  have  a  layered  structure,  an  internal  matric  crack 

» 

cannot  always  undergo  monotonic  propagation;  the  crack  may  either  be  arrested 
or  blunted  at  some  layer  Interface  boundaries.  This  crack-arresting  mech¬ 
anism  permits  an  Increase  of  the  applied  load  without  causing  catastrophic 
failure  of  the  laminate.  Consequently,  other  flaws  in  the  laminate  may 
become  critical  and  similar  matrix  cracks  are  formed  as  the  applied  load 
is  Increased. 

• 

Thus,  after  a  certain  loading  process,  be  it  monotonic  or  cyclic 
stressing,  a  distribution  of  matrix  cracks  within  the  laminate  member  is 
obtained.  • 

Although  the  inherent  material  microflaws  may  exist  randomly  in  the 
laminate,  the  general  characteristics  of  the  macroscopic  matrix  cracks 

are  found  to  be  highly  dependent  on  the  nature  of  loading  as  well  as  the  * 

lamination  geometry,  such  as  laminate  stacking  sequence,  lamina  thickness. 

laminate  shape  features,  etc.  [2].  Generally,  for  a  given  laminate  under  a 

given  class  of  loading,  the  matrix  cracking  pattern  is  mechanically  repro-  * 

duclble,  though  it  may  exhibit  a  certain  statistical  variation  [3]. 
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This  fact  suggests  that  the  formation  mechanisms  of  a  certain  mode 
of  sub-laminate  cracks  is,  at  least  partially,  mechanically  deterministic, 
while  the  precursors  of  the  cracks  (material  microflaws)  remain  probabil¬ 
istic  . 

This  observation  has  led  the  development  of  the  desired  cumulative 
damage  model  to  a  methodology  which  combines  both  a  probabilistic  and  a 
mechanistic  approach. 

The  actual  procedures  followed  for  the  model  construction  in  this 
research  program  is  illustrated  by  the  schematic  flow  diagram  shov«m  in 
Figure  1.1.  The  various  developmental  stages  in  this  flow  diagram  are  dis¬ 
played  in  a  two-dimensional  grid-work  as  shown.  From  the  top  to  the  bottom 
row  there  are  three  parallel  tasks;  namely: 

(a)  Experimental  investigation  on  laminates  in  which  a  certain  sub¬ 
laminate  cracking  mode  is  developed  under  external  loading. 

The  entire  load  and/or  time  dependent  crack  growth  process  is 
monitored  and  documented.  The  physical  parameters  that  influ¬ 
ence  the  crack  development  mechanisms  are  identified  and 
analyzed . 

(b)  Analytical  investigation  in  which  the  development  of  the  specific 
sub-laminate  cracking  process  is  simulated  by  a  stochastic  pro¬ 
cedure  incorporating  the  mechanics  of  fracture  and  the  concept  of 
random  material  flaw  distributions. 

(c)  A  correlation  between  the  experiment  and  the  analytical  simula¬ 
tion  whereby  the  desired  damage  cumulative  model  is  constructed. 


In  each  of  the  three  tasks  outlined  above,  there  are  four  investi¬ 
gative  stages  shown  in  Figure  1.1;  from  left  to  right  they  are: 

(a)  A  study  of  the  damage  accumulation  processes  under  static 
loading.  In  this  stage,  the  major  modes  of  sub-laminate  cracking 
are  identified.  In  each  case,  the  predominant  influencing 
physical  parameters  are  taken  into  account  in  the  simulation 
model  and  a  load -damage  growth  relation  is  then  established. 

(b)  A  study  of  the  damage  accumulation  processes  under  constant  ampli¬ 
tude  cyclic  loading.  In  this  stage,  a  crack  growth  law  is  de¬ 
veloped  and  a  load-cycle-damage  relation  is  established  for  each 
of  the  major  failure  modes. 

(c)  From  the  results  of  stages  (a)  and  (b)  above  a  relationship 
between  the  load  (o)  and  the  fatigue  cycle  (N)  is  established 
for  a  given  laminate  damage  state. 

(d)  Using  the  results  of  (c) ,  a  model  for  the  cumulative  damage 
process  in  laminates  subjected  to  variable  amplitude  fatigue 
loading  is  then  constructed  based  on  the  concept  of  the  constant 
damage  states. 

In  Chapter  II  of  this  report  we  shall  present  the  details  of  the  ana¬ 
lytical  developments  as  highlighted  above.  The  predictive  models  are  then 
applied  to  the  specific  cases  for  which  experimental  results  have  been 
generated. 

Chapter  III  describes  the  scope  of  the  experiment.  Most  of  the  relevant 
results  are  either  tabulated  and/or  graphically  displayed  in  the  Appendices. 


A  brief  discussion  on  the  current  status  of  the  research  program 
and  a  summary  of  the  continued  work  in  Phase  III  are  highlighted  in 
Chapter  IV. 

• 

» 

» 
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PAHACE  MODEL  CONSTRUCTION 


II.  ANALYTICAL  DEVELOPMENT 


2.1.  Major  Modes  of  Sub-Laminate  Cracking 

Experiments  using  graphite-epoxy  composite  laminates  have  shown  that 
transverse  cracks,  fiber  splittings  and  free-edge  induced  delamination 
are  among  the  most  frequently  encountered  sub-laminate  cracking  modes. 
Furthermore,  these  cracking  modes  can  occur  independently  or  concurrently 
depending  on  the  individual  case  conditions.  Sometimes,  one  mode  can 
Interact  with  another  resulting  in  a  complicated  damage  state  in  the  lami¬ 
nate. 

For  purpose  of  an  analytical  model  construction,  it  is  necessary  to 
first  understand  the  exact  mechanisms  of  each  individual  cracking  mode 
and  then  the  mechanisms  of  two  or  more  Interacting  modes.  In  this  research 
program,  it  has  been  planned  that  the  following  cracking  mechanisms  be 
investigated : 

(a)  Initiation  and  growth  of  multiple  transverse  cracks  in  lami¬ 
nates  containing  a  90°  layer, 

(b)  Initiation  and  growth  mechanisms  of  a  predominantly  free-edge 
Induced  interlaminar  cracking  (delamination) , 

(c)  Interaction  between  multiple  transverse  cracking  and  free-edge 
delamination,  and 

(d)  Interaction  between  multiple  transverse  cracking  and  longitudinal 
(0° -layer)  fiber  splitting  mechanisms. 

In  all  cases  both  the  static  growth  mechanisms  and  the  fatigue  growth 
mechanisms  are  to  be  studied.  However,  since  Items  (a)  and  (b)  above 
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require  a  two-dimensional  stress  analysis,  while  items  (c)  and  (d)  require 
a  three-dimensional  stress  analysis,  the  present  research  effort  follows 
the  same  order  as  listed  above. 


In  this  chapter,  a  unified  methodology  is  to  be  developed  which  is 
based  on  the  fundamental  concept  of  the  classical  fracture  mechanics  and 
the  concept  of  effective  material  flaws. 

2.2.  Application  of  Fracture  Mechanics 

To  facilitate  our  subsequent  model  development,  let  us  first  intro¬ 
duce  briefly  the  basic  concept  of  the  classical  fracture  mechanics. 
Consider  the  well-known  Griffith  problem  which  pertained  to  an  elastic 
homogeneous  and  isotropic  plate  uniformly  stretched  in  one  direction  by  o. 
The  plate  has  a  through  crack  of  size  2a  orientated  normal  to  the  direction 
of  a.  The  length  of  the  crack  is  assumed  small  but  finite.  Griffith  pos¬ 
tulated  that:  at  the  instance  of  crack  extension,  a  loss  of  the  stored 
elastic  strain  energy  near  the  crack-tip  region  is  incurred.  This  energy 
is  assumed  to  be  convered  into  the  surface  energy  of  the  crack  extension. 
Thus,  a  balance  of  energy  at  the  instance  of  crack  extension  leads  to  the 
criterion 

m 

3a  3a  ^ 

Equation  (1)  defines  the  general  condition  under  which  the  existing 
crack  begins  to  propagate  in  self-similar  mode. 

For  the  plate  considered  by  Griffith,  the  critical  stress  at  the  in¬ 
stance  of  crack  extension  can  be  derived  from  (1),  giving 


(2) 


a 

c 


( 


2Ey 

wa 


1/2 

) 


where  y  is  the  free  surface  energy  density  (per  unit  area)  of  the  material. 


When  applied  to  engineering  problems,  the  Griffith  theory  is  often 
modified  for  practical  considerations.  For  instance,  in  Equation  (2) 
becomes  unbounded  as  a  -»■  0.  Clearly,  no  real  material  can  sustain  an  in¬ 
finite  stress.  The  limitation  of  Equation  (2)  can  be  conveniently  circum¬ 
vented  by  invoking  the  assumption  of  material  flaws.  That  is  to  assume 
that  the  material  possesses  naturally  some  characteristic  distribution  of 
microflaws.  It  is  then  assumed  that  the  worst  flaw,  having  a  finite 
empirical  size  of  2a^,  acts  like  a  real  crack.  This  particular  flaw  de¬ 
termines  a  finite  critical  stress  according  to  Equation  (2) . 

Of  course,  flaws  do  exist  in  real  materials,  especially  in  fibrous 
composite  systems.  But,  the  real  physical  identity  of  the  flaws  cannot 
be  practically  Included  in  the  fracture  analysis;  only  their  aggregate 
effects  could  perhaps  be  considered.  Thus,  to  apply  the  concept  of  the 
classical  fracture  mechanics  to  problems  involving  cracks  in  composites, 
a  new  concept  of  "effective  material  flaw  distribution"  will  be  introduced. 
Detailed  discussion  of  this  concept  is  presented  in  the  next  section. 


Another  practical  consideration  is  related  to  the  definition  of  y, 
the  free  surface  energy  density  of  the  material.  For  a  crack  in  a  brittle 
material  such  as  glass,  y  could  be  considered  as  defined  above.  For  most 
other  engineering  materials,  however,  crack  extension  is  found  to  associate 
a  certain  degree  of  inelastic  deformation  near  the  crack-tip  region.  For 
this  reason  a  crack  extension  will  consume  more  than  just  the  elastic 
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strain  energy  released;  It  may  also  consume  some  amount  of  the  inelastic 
deformation  energy.  Furthermore,  the  crack  extension  path  or  the  crack 
surface,  may  show  a  certain  degree  of  ruggedness  depending  on  the  hetero¬ 
geneity  of  the  material  velwed  at  a  microscale.  Early  studies  on  struc¬ 
tural  metals  have  regarded  the  quantity  y  as  the  total  irreversible  work 
required  to  create  a  unit  crack  surface  area.  Or,  it  can  be  interpreted 
as  the  energy  dissipated  in  the  crack-tip  region  during  crack  extension. 
Then,  y  will  depend  on  the  inelasticity  and/or  the  microscopic  hetero¬ 
geneity  of  the  material  locally.  It  is  regarded  also  as  an  intrinsic  prop 
erty  of  the  material  which  can  be  defined  only  empirically. 

These  considerations  are  of  fundamental  importance  when  a  crack-like 
failure  in  composite  laminate  is  to  be  modeled  using  fracture  mechanics 
method . 

In  common  practice  the  right-hand  side  of  Equation  (1)  is  replaced 
by  =  (2y) ,  known  as  the  critical  energy  release  rate  of  the  material. 
The  left-hand  side  is  a  function  of  the  applied  load,  the  geometry  of  the 
body,  and  the  size  of  the  crack.  Thus,  for  a  crack  undergoing  self- 
similar  extension,  the  Griffith  criterion  is  expressed  as 

G (o , a)  =  G  •  ( ^ 

c 

Accordingly,  the  development  of  the  fracture  model  for  a  known  sub- 

laminate  crack  rests  upon  the  calculation  of  G(a,a)  and  the  physical 

measurement  of  G  . 

c 

As  has  been  mentioned  earlier,  the  quantity  a  in  G(a,a)  represents 
the  effective  flaw  size  which  is  actually  random  in  nature.  Hence,  a  new 
concept  must  be  developed  which  includes  this  Important  feature. 
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The  Concept  of  Effective  Material  Flaw  Distribution 

In  stress  analysis  of  laminates,  the  basic  approach  has  been  the 
concept  of  ply-elastlclty  [4].  This  concept  pertains  to  laminates  that 
are  made  of  unidirectional  composite  plies  and  are  stacked  together 
with  each  ply  orientated  in  a  certain  designed  direction.  For  purpose 
of  laminate  analysis,  the  unidirectional  ply  is  regarded  as  the  building 
block  of  the  laminate.  The  approach  is  to  characterize  the  effective 
properties  of  the  building  block  so  as  to  predict  the  properties  of  the 
laminate.  The  fundamental  assumption  in  ply-elasticity  is  that  the  uni- 
dlrectlonaly  ply,  which  is  microscopically  heterogeneous,  can  be  replaced 
by  an  "effective"  homogeneous  body  without  altering  the  global  properties 
of  the  laminate.  The  assumption  thus  enables  the  representation  of  the 
unidirectional  ply  by  two  seta  of  homogeneous  constants.  The  first  is 
associated  with  the  deformation  properties,  e.g.,  via  the  generalized 
Hookes  Law: 


(4) 


where  the  stiffness  constants,  ,  are  determined  as  averaged  quantities 
over  a  solid  volume  of  macroscopic  dimension. 

The  second  set  of  material  constants  is  associated  with  the  failure 
properties,  e.g.,  via  the  strength  criterion. 


(5) 


where  and  F^  are  strength  constants  determined  upon  the  first  failure 
of  the  ply  under  the  homogeneous  stress  condition  o^. 

However,  it  is  generally  accepted  that  failure  o^  structural  material 
is  a  local  extremum  property  [5].  The  tensile  strength  in  the  transverse 
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direction  of  the  unidirectional  ply,  for  Instance,  is  determined  by  some 
local  extrema.  Local  extrema,  of  course,  included  the  inherent  material 
microflaws  which  form  a  distribution  both  in  size  and  in  location.  It 


is  clear  then,  the  ply  strength  constants  and  are  determined  essen¬ 
tially  by  the  "dominating"  local  flaw  in  the  unidirectional  ply*  The  pre¬ 
sence  of  other  (smaller)  flaws  cannot  be  reflected  by  the  constants  F^^ 


and  F 


ij' 


In  order  to  model  multiple  failures  in  laminate  such  as  described 
earlier,  it  is  necessary  to  Include  the  existence  of  the  microflaw  dis¬ 
tribution  in  the  ply.  Hence,  it  will  be  postulated  here  that  multiple 
failures  in  a  lamina  are  governed  by  the  Initiation  and  propagation  of  a 
distribution  of  "effective"  flaws.  These  hypothetical  flaws  act  indi¬ 
vidually  like  cracks  and  are  capable  of  propagating  into  matrix  cracks 
during  a  certain  course  of  loading.  In  this  manner,  failures  in  the 
lamina  will  be  governed  by  a  criterion  of  crack  propagation  rather  than 
by  a  homogeneous  failure  law  such  as  Equation  (5). 


Specifics  of  the  effective  flaw  distribution  will  be  covered  when 
specific  problems  are  discussed.  At  this  point,  let  us  just  regard  the 
quantity  a  in  Equation  (3)  as  a  given  random  variable  for  the  material 
system. 

The  Calculation  of  G(q,a) 

In  the  theory  of  the  LEFM,  the  singular  stress  field  near  a  crack  tip 
in  a  homogeneous,  isotropic  elastic  body  has  been  represented  by  some 
analytical  functions  in  the  theory  of  complex  variables  [6].  These  near¬ 
field  stresses  can  result  in  three  particular  modes  of  the  crack  opening. 


These  are  known  as  the  opening  mode  (I),  the  sliding  mode  (II),  and  the 
anti-plane  shearing  mode  (III) .  For  each  mode,  the  stresses  are  expressed 
in  terms  of  the  associated  stress  intensity  factor  K  [7];  and  consequently, 
the  associated  G(o,a)  is  computed  in  terms  of  K.  Since  the  relation 
between  K  and  G  is  one-to-one.  Equation  (3)  reduces  to  the  form 

K(o,a)  =  (6) 

Similar  relations  between  K  and  G  for  orthotropic  media  having  a  crack 
orientated  along  one  of  the  major  axes  can  also  be  obtained  [8].  But  the 
analytical  solutions  for  the  singular  stress  field  often  require  tedious 
mathematical  derivations. 

Direct  solution  methods  for  G(a,a)  have  been  available,  among  them  are 
the  well-known  J-integral  method  [9,  10)  and  the  virtual  crack  closure 
technique  by  Irwin  [11). 

Irwin  [11)  observed  that  the  elastic  strain  energy  released  during  a 
virtual  crack  extension,  Aa,  is  equal  to  the  work  done  in  closing  it  again. 
But  the  Inverse  problem  provides  easily  the  solution  for  the  surface 
transactions,  5,  over  a.  The  crack-tip  energy  release  rate  is  then  rep¬ 
resented  by 

<=  ■  mC  <■“  ■ 

where  Au  is  the  crack  opening  displacement  vector  over  Aa. 

If  the  crack  extension  involves  all  three  modes  (I,  II,  III),  the 
vector  product  in  Equation  (7)  will  give  a  sum  of  three  scalars  associated 
respectively  with  G^,  and 
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The  virtual  crack-closure  representation  is  particularly  adaptive 
to  machine  computations.  Ryblcki  and  Kanninen  [12]  suggested  a  two- 
dimensional  finite  element  technique  to  evaluate  G  for  a  line  crack  in  a 
plane.  The  crack-tip  stress  vector,  o,  and  the  displacement  vector,  &u, 
in  Equation  (7)  are  approximated  by  the  nodal  forces  and  displacements, 
respectively,  in  a  finite  element  representation. 

Wang  and  Crossman  applied  this  technique  and  developed  a  generalized 
plane  strain  finite  element  routine  [13,  14],  which  can  similate  a  line 
crack  propagation  In  the  two-dimensional  cross-section  of  a  laminate. 

Since  a  general  laminate  under  load  may  suffer  cross-sectional  warping, 
the  routine  actually  computes  three-dimensional  stresses  and  displacements 
and  hence  also  the  three  modes  of  the  energy  release  rate  G^,  and 

If  a  laminate  is  subjected  to  the  far-field  stress  and  a  sub- 
laminate  crack  is  to  be  simulated,  the  crack-tip  energy  release  rate  can 
be  expressed  in  the  general  form 

G  =  C  (a/t)  •  t  .  (o  /E  (8) 

e  e  o  o 

where  E  is  laminate  stiffness  in  o  direction  and  t  is  the  linear  scale 
o  o 

between  the  actual  model  and  the  finite  element  model.  C  is  a  function 

e 

dependent  only  on  the  crack  size  a  which  is  routinely  generated  for  a 
known  sub-laminate  cracking  mode  in  a  given  laminate,  see  [13]. 

Similarly,  if  the  laminate  is  subjected  to  a  uniform  temperature  drop 
At  and  if  a  thermally  induced  crack  opening  is  resulted  then  the  associ¬ 
ated  energy  release  rate  at  the  crack  tip  is  expressed  by 

G.|,  =  C^(a/t)  •  t  •  (AT)^  (9) 


14 


I 


Generally,  the  laminate  is  prestressed  by  a  AT  due  to  curing,  so 
a  combined  effect  is  resulted  when  is  applied;  the  combined  energy 


release  rate  is: 


G  =  [(C  e  +  (C  AT]^ 

e  o  1 


where  e  is  the  far-field  strain  (=  o  /E  ) . 
o  o  o 

It  is  seen  that  the  finite  element  technique  is  extremely  versatile, 
and  can  be  efficiently  executed  for  simulating  complicated  sub-laminate 
cracking  growth  such  as  delamlnatlon. 

In  this  research  program  the  finite  element  routine  developed  by 
Wang  and  Crossman  [13,  14]  was  applied  for  the  energy  release  rate  cal¬ 
culations.  Some  details  of  this  routine  have  also  been  discussed  in  the 
Phase  I  Final  Report  [15]. 


The  Evaluation  of  G 


When  fracture  occurs  in  the  material,  the  energy  released  in  the  pro¬ 
cess  is  expected  to  depend  in  part  on  the  microscopic  morphology  of  the 
fracture  surface.  Fibrous  composites  are  known  to  possess  very  complex 
fracture  surface  details  [16].  The  observed  delamination  surfaces  in 
graphite-epoxy  composites,  for  instance,  reveal  a  considerable  raggedness 
because  the  crack  must  pass  around  the  reinforcing  fibers.  But,  any  such 
crack  surface  details  will  be  "smoothed  out"  by  the  assumption  used  in 
ply-elasticity  and  their  effects  can  be  reflected  only  in  the  measured 
quantity  of  the  fracture  energy  release  rate  G^. 

For  this  reason,  G^  measured  for  some  matrix-predominant  fractures 
in  epoxy-based  composites  has  been  found  to  depend  on  the  direction  of 
fracture  propagation.  Cullen  [16]  and  Williams  [17]  considered  two 
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different  cases  of  delamination  in  a  graphite-epoxy  composite  system 

as  illustrated  in  Figure  2.1.  The  first  case  is  a  0°/0°  delamination  I 

in  which  the  crack  path  is  in  the  fiber  direction;  while  the  second  is 
a  90°/90°  delamlnatlon  where  the  crack  path  is  transverse  to  the  fiber 

direction.  Both  were  subjected  to  mode-I  cracking  condition.  They  found  I 

that  the  microscopic  morphology  of  the  90®/90°  delamination  surface  ex¬ 
hibited  considerably  more  raggedness  than  the  O^/O®  delamination  surface. 

This  resulted  in  a  marked  difference  for  the  measured  G  .  Note  that  I 

c 

these  two  fracture  events  occur  essentially  in  the  same  interface  when 

viewed  macroscopically .  Yet,  the  respective  values  can  differ  greatly 

depending  on  the  direction  of  the  crack  propagation.  I 

Clearly,  when  the  Griffith  formula  Equation  (3)  is  applied  for  crack 
extension  in  composites,  the  term  G^  requires  a  precise  qualification. 

Similarly,  when  a  test  method  is  devised  to  measure  G^,  the  dimensional  and  ^ 

directional  characteristics  of  the  measured  data  must  also  be  considered. 


Mode-I  Interlaminar  G.,  .  A  commonly  used  test  method  to  determine  inter- 

laminar  G^,^  is  the  splitting  cantiliver  beam.  Cullen  [16]  and  Wilkins  [18] 

have  used  this  method  to  determine  the  interlaminar  G,.  when  the  crack  is 

Ic 

propagating  in  the  direction  of  the  fibers  (0°/0°  delamination).  For  the 

graphite-epoxy  systems  used  in  their  study,  (AS-3501-02) ,  they  found  that 

2 

Gj^  at  room  temperature  is  about  0.85  Ib/in,  or  130  j/m  .  This  value  is 
about  twice  the  G^^  measured  for  the  pure  epoxy  resin.  Williams  [17]  used 
a  compact  specimen  which  simulates  roughly  a  90°/90°  delamination.  He 
found,  for  the  same  material  system,  a  G^^  value  of  1.3  Ib/in,  or  about 
three  times  that  of  the  pure  resin.  Williams  explained  that  the  fracture 
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surface  In  his  specimens  showed  fiber  breakage  as  well  as  fractured 


epoxy  debris;  this  had  resulted  in  a  higher  value  for 

In  another  paper  by  Wilkins,  et  al.  [19],  it  is  reported  that 

in  delamination  of  0°/90°  interface  is  also  found  higher  than  that  of 

0°/0°  delamination.  These  findings  reaffirm  the  directional  dependent 

nature  of  G,.  . 

Ic 

Mixed-mode  interlaminar  G^_  .  The  splitting  cantiliver  beam  method 

- (I,  II)c 

has  also  been  used  in  mixed-mode  crack  growth  experiments.  In  this  case, 

it  is  necessary  to  apply  different  loads  at  the  top  and  the  bottom  parts 

of  the  split  so  as  to  create  both  an  opening  (1)  and  a  sliding  (II) 

action.  Vanderkley  [20]  and  Wilkins  [18]  conducted  tests  on  the  same 

graphite-epoxy  system  (used  for  their  mode-I  tests)  and  found  the  total 

energy  release  rate  G._  =  G_  +  G,,  that  exhibited  a  strong  depen- 

(I,  II)c  Ic  IIc 

dence  on  the  G^^/G^  ratio. 

This  phenomenon  (G^j/Gj^  ratio  dependence)  is  not  uncommon  in  a  mixed¬ 
mode  fracture.  Similar  observations  were  reported  for  brittle  metals  as 
well  as  pure  epoxy  resins  [21].  It  is,  perhaps,  more  pronounced  in  fibrous 
composites.  Generally,  the  increased  fracture  resistance  in  mixed-mode 
cracks  is  attributed  to  possible  matrix  yielding  under  shear  as  well  as 
crack-closure  friction  due  to  the  sliding  action.  In  any  event,  this 
question  requires  further  research. 

In  a  test  on  a  series  of  double-notched  off-axis  directional  graphite- 

epoxy  (AS-3501-06)  laminates,  Wang,  et  al.  [22]  measured  the  mixed-mode 

G(i  a  continuous  function  of  Gj^/G^;  their  results  are  shown  in 

Figure  2.2.  It  is  seen  that  G.^  is  monotonically  increasing  with 

(I,II)c 
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Gjj/Gj,  possibly  reaching  an  asymtotic  value. 

Having  established  a  reasonable  understanding  about  the  quantities 
on  both  sides  of  the  fracture  equation.  Equation  (3),  we  are  now  ready 
to  apply  this  equation  to  some  specific  types  of  sub-laminate  cracking. 
The  remainder  of  this  chapter  is  devoted  entirely  for  this  purpose. 


I 

2.3.  Mechanics  of  Multiple  Transverse  Cracks  -  Static  Loading. 

Transverse  cracks  are  found  in  epoxy-based  laminates  usually  at  a  * 

relatively  low  loading  level.  When  viewed  at  the  macroscopic  scale,  the 
cracking  action  is  simply  a  sudden  separation  of  fibers  by  breaking  the 
epoxy  bond.  To  illustrate,  consider  as  an  example  a  [0/90]^  type  lami-  ^ 

nate  such  as  shown  in  Figure  2.3.  Under  the  far-field  tensile  stress 
the  90®-layer  may  suffer  multiple  transverse  cracks.  Generally,  the 
sequence  of  events  is  as  follows:  a  single  crack  forms  first  when  the  * 

applied  load  reaches  a  certain  critical  value,  which  defines  the  "onset" 
of  the  events;  as  the  applied  load  increases,  more  similar  cracks  are 
formed.  If  there  is  no  other  failure  mode  setting  in  at  high  load,  the  ^ 

number  of  transverse  cracks  will  continue  increasing,  until  it  reaches 
a  saturation  density. 

Figures  2 . 4  and  2.5  show  the  experimentally  obtained  crack  density 
(number  of  transverse  cracks  per  unit  length  of  specimen)  versus  the 
applied  far-field  stress  plots  for,  respectively,  [^2^^®2^s  ^^2^^^3^s 

laminates  (these  results  were  reported  earlier  in  Reference  [15]).  It  is 
seen  that  the  thickness  of  the  90°-layer  in  these  two  laminates  has  a 
profound  influence  on  the  characteristics  of  the  crack  formation  and 

I 

growth  behaviors.  Specifically,  the  laminate  with  thinner  90°-layer  re¬ 
quires  a  higher  load  for  onset  of  transverse  cracking;  it  is  also  capable 
of  yielding  a  higher  crack  density  than  the  laminate  with  the  thicker 

l 

90°-layer . 

The  90°-layer  thickness  effect  on  transverse  cracking  behavior  in 
[0/90]^  laminates  was  first  documented  experimentally  by  Baser,  et  al. 

[23].  The  effect  was  attributed  to  the  constraining  actions  of  the 
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adjacent  0°-layers;  they  observed  that  a  transverse  crack  can  be  no 
larger  than  the  90°-layer  thickness  and  the  energy  release  rate  which 
drives  the  crack  is  also  limited  by  the  same  factor. 

The  90°-layer  thickness  effect  on  crack  density  has  been  explained 
by  the  existence  of  a  "shear-lag"  zone  at  the  transverse  crack  root  [23  - 
25].  That  is,  an  interlaminar  shear  stress  zone  is  developed  near  the 
0/90  Interface  where  a  transverse  crack  terminates.  This  shear  stress 
is  singular  at  the  crack  root  but  decays  exponentially  a  distance  away. 
Inversely,  the  in-situ  tensile  stress  in  the  90°-layer  is  nil  on  the 
transverse  crack  surface  but  it  regains  its  far-field  value  outside  the 
shear-lag  zone.  Ideally,  any  two  adjacent  cracks  should  be  spaced  by 
the  shear-lag  distance  and  the  shear-lag  distance  is  proportional  to  the 
90°-layer  thickness,  hence  the  observed  90®-layer  thickness  effect  on 
crack  density. 

Although  the  shear- lag  argument  is  Ideally  correct,  an  evenly  spaced 
transverse  cracks  pattern  does  not  occur  in  practice;  and  the  transverse 
cracking  behavior  is  Influenced  strongly  also  by  the  Inherent  flaw  dis¬ 
tribution  in  the  90°-layer. 

The  Energy  Model  for  Transverse  Cracks 

For  purpose  of  clarity,  consider  a  [0/90]^  type  laminate  such  as 
shown  in  Figure  2.3.  Assume  that,  in  the  90'’-layer,  the  material  has  a 
random  distribution  of  microflaws.  The  gross  effects  of  the  microflaws 
at  the  macroscopic  scale  are  represented  by  a  liypothet i ea 1  distribution 
of  "effective"  flaws  which  cannot  by  physically  identified  at  the  macro- 
>  scale.  But,  under  stress,  these  etfective  flaws  are  capable 
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of  propagating  into  transverse  cracks,  which  are  physically  real.  The 
individual  size  of  the  effective  flaws  is  denoted  by  2a  and  any  two  ad¬ 
jacent  flaws  are  spaced  by  a  distance  S,  see  Figure  2.6.  Let  us  then 
assume  that  for  a  given  laminate  specimen  of  a  unit  length,  a  character¬ 
istic  effective  flaw  distribution  can  be  found  and  it  is  represented  by 
the  probability  density  function  f(a)  and  the  probability  density  func¬ 
tion  f(S).  For  sake  of  simplicity  and  of  no  evidence  to  suggest  otherwise, 
the  two  functions  are  assumed  to  take  the  form  of  a  normal  distribution, 


f  (a) 

J. 

exp  [-  (a  -  p  )^/2o 
a  a 

(11) 

f(S) 

1 

s/27 

exp  (-  (S  -  p  )^/2a 
s  s 

(12) 

where  h  and  <y  are  the  mean  and  the  standard  deviation  of  the  respective 
distribution  functions. 

Among  all  the  "effective"  flaws,  let  the  size  of  the  "worst"  one  be 
denoted  by  2a  .  Then,  this  flaw  determines  the  first  failure  of  the  90°- 
layer.  The  propagation  of  the  first  flaw  under  the  applied  laminate  load 
thus  forms  the  "onset"  of  the  transverse  cracking  process.  With  increased 
applied  load,  flaws  of  lesser  size  continue  to  form  transverse  cracks 
allowing  multiple  failures  of  the  SO^-layer. 

In  the  representations  (11)  and  (12),  the  parameters  u  ,  a  ,  p  ,  and 
Og  are  not  ,  strictly  speaking,  properties  of  a  unidirectional  ply. 

Rather,  it  is  more  appropriate  to  regard  these  as  properties  of  the  90°- 
layer  in  the  laminate  under  consideration.  For  example,  the  size  of  the 
"worst"  flaw  2a^  must  be  smaller  than  the  thickness  of  the  90°-layer,  2h, 
see  Figure  2.6.  Thus,  for  a  normal  distribution  such  as  expressed  in  (11), 
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let  a  be  the  99  percentile  of  f(a).  Then,  the  following  must  be  valid 
o 

[26]: 

a  -  p  +  30  <  h  (13) 

o  a  a 


Note  that  h  can  be  made  arbitrarily  small  or  large  depending  on  the  lami¬ 
nation  design. 


Within  the  framework  of  the  foregoing  discussions,  the  "effective" 
flaw  concept  becomes  readily  amenable  to  methods  of  classical  frac¬ 
ture  mechanics.  To  illustrate,  consider  the  laminate  shown  in  Figure  2.6. 
A  distribution  of  "effective"  flaws  is  assumed  to  exist  along  the  length 
of  the  90°-layer.  Let  the  distribution  be  represented  by  the  probability 
functions  f(a)  and  f(S)  as  expressed  in  (11)  and  (12)  and  the  parameters 
p  and  0  in  these  functions  are  assumed  known.  Then,  the  "worst"  flaw 
whose  size  is  2a^  will  propagate  and  form  the  first  transverse  crack  when 
the  available  energy  release  rate,  G(a^,  ’  equals  the  luaterial's  criti¬ 

cal  energy  release  rate  G^: 


G(0,a)  —  G. 

X  o  c 


The  stability  of  the  propagation  is  stable  if 


(14) 


G(o  ,  a  +  Aa)  <  G 
X  o  c 


(15) 


unstable  if 


G(o  ,  a  +  A  )  >  G  . 
X  o  a  c 


(16) 


Under  the  laminate  tensile  stress,  o  ,  G  is  a  function  of  the  effec- 

X 

tive  flaw  size  a.  The  general  shape  of  this  function  is  Illustrated  in 
Figure  2.7.  Calculation  of  G  for  a  specific  laminate  may  be  performed 
by  a  numerical  procedure  based  on  ply  elasticity  (see  details  in  Ref.  [13]) 
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The  critical  laminate  load  at  the  onset  of  the  first  crack  is  deter¬ 
mined  from  the  fracture  criterion  (14),  given  the  value  of  a^. 

The  second  crack  will  form  from  another  flaw  in  the  flaw  distribution 
f(a)  when  the  same  energy  condition  surrounding  that  flaw  is  reached. 

But,  this  flaw  not  only  is  smaller  in  size  but  also  is  located  at  a  random 
distance  from  the  first  crack.  Assuming  the  far-field  laminate  load 
remains  unchanged,  the  stresses  near  any  flaw  which  is  located  in  the  close 
vicinity  of  the  crack  will  be  reduced.  The  available  energy  release  rate 
at  the  tip  of  the  flaw  depends  on  its  relative  distance  S  from  the  crack, 
see  Figure  2.8.  This  quantity  may  be  expressed  in  the  following  form; 

G(a  ,  a)  =  R(S)  •  G  (5  ,  a)  (17) 

X  ox 

where  G^  is  the  energy  release  rate  of  the  flaw  of  size  2a  without  the  pre¬ 
sence  of  a  neighboring  crack,  and  R(S)  is  the  energy  release  rate  retention 
factor  when  the  presence  of  the  neighboring  crack  is  considered.  Depending 
on  the  value  of  S,  R(S)  may  take  a  value  between  0  and  1,  Figure  2.8. 

Hence,  given  the  values  of  a  and  S,  the  laminate  load  corresponding 
to  the  formation  of  the  second  crack  can  be  determined  from  the  energy 
criterion: 

R(S)  ,  a)  =  G  (18) 

Ox  c 

Generally,  ,  determined  from  (18),  is  slightly  higher  than  that  de¬ 
termined  for  the  first  crack. 

For  a  flaw  located  between  two  cracks,  the  energy  release  rate  is 
given  by: 

G(a  ,  a)  =  R(S  ).G  (a  ,  a).R(S  )  (19) 

X  J6  O  X  X* 

where  and  S^,  are  respectively,  distances  to  the  left  and  to  the  right 
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cracks  from  the  considered  flaw. 

To  determine  the  entire  load-sequence  of  the  multiple  crack  formation 
process,  a  stochastic  simulation  procedure  will  be  followed.  The  procedure 
essentially  simulates  the  cracking  process  as  it  would  occur  naturally. 

Each  simulation  represents  a  test  case  of  a  sample;  and  repeated  simulations 
represent  tests  on  a  group  of  replicates. 

A  Stochastic  Simulation  Procedure. 

The  simulation  is  basically  a  numerical  Monte-Carlo  random  search 
routine  [27].  A  flow  diagram  of  the  computer  routine  is  illustrated  in 
Figure  2.9.  The  routine  combines  both  the  fracture  mechanics  and  the 
probabilistic  representation  of  the  material  flaw  distribution.  Specifi¬ 
cally,  fracture  mechanics  provides  a  criterion  for  the  flaws  to  become 
cracks;  while  the  probabilistic  flaw  distribution  yields  a  continuous 
load-sequence  for  the  cracks  formation  process. 

The  flow  diagram  in  Figure  2.9  shows  that  input  information 
is  generated  for  the  specific  laminate  which  is  to  undergo  crack  simula¬ 
tion.  First,  the  energy  release  rate  a)  is  generated  for  a  flaw 

of  2a  without  the  presence  of  a  neighboring  crack.  Both  and  a  are 
treated  as  arbitrary  parameters.  This  is  done  by  expressing 
the  form  [2] : 

G  (o  ,  a)  =  •  e  +  /cl  •  AT]^  .  t  (20) 

ox  e  o  T 

where  e  is  the  laminate  strain  (=  o  /E  ) ,  AT  denotes  the  uniform  tempera- 

O  XX 

ture  drop  due  to  laminate  curing,  t  is  a  characteristic  length  parameter, 
and  and  are  coefficient  functions  that  depend  only  on  a.  Thus,  a 
separate  calculation  of  and  is  performed  before  the  simulation. 


Secondly,  the  energy  release  rate  retention  factor ,  R(S),  is  generated 
for  the  considered  laminate  by  computing  G(0^,  a)  for  a  flaw  of  size  2a 
under  the  Influence  of  a  transverse  crack  located  at  a  distance,  S,  away, 
(see  Figure  2.8). 

Another  input  to  the  simulation  routine  is  the  knowledge  of  the  lami¬ 
nate  ultimate  strength  a^.  This  is  used  to  merely  terminate  the  simula¬ 
tion  once  the  applied  load,  a  ,  reaches  a  . 

With  the  above  input  information,  the  simulation  starts  with  gener¬ 
ating  a  set  of  n  uniformly  random  values  in  the  Interval  (0,1).  Then,  by 
equating  the  cumulative  function  F(a)  of  the  assumed  flaw  size  distribution 
function  f (a)  to  each  of  the  random  values,  a  set  of  n  values  for  a  is  com¬ 
puted,  say,  {a^},  i  =  l,n.  Here,  one  needs  to  assign  a  priori,  appropriate 

values  for  p  and  a  ,  in  order  to  computer  {a  }.  These  parameters  will  be 
3d  X 

correlated  later  with  experiments. 

Similarly,  a  set  of  (n-1)  values  for  S  is  also  computed,  {S^},  i  =  1, 
(n-1) .  is  then  assigned  as  the  spacing  between  the  ith  and  the  (i+l)th 

flaws. 

The  computer  search  begins  with  the  identification  of  the  flaws 
most  likely  to  become  transverse  cracks.  The  first  to  occur,  clearly,  is 
the  largest  flaw  in  the  ta^}  set.  The  corresponding  applied  load  for  the 
onset  of  the  first  crack  in  the  process  is  determined  by  the  fracture 
criterion  (14) . 

After  the  first  crack,  the  routine  searches  for  the  next  flaw  most 
likely  to  become  a  crack  by  comparing  the  G(o^,  a)  values  of  all  the 
other  flaws  with  the  presence  of  the  first  transverse  crack.  Equation  (18) 


is  used  to  determine  the  applied  load  associated  with  the  next  crack.  Simi 


larly,  for  the  rest  of  the  cracks,  the  corresponding  loads  are  determined 

using  (19).  The  simulation  goes  on  until  the  applied  load,  a^,  reaches  the 

laminate  strength,  o  . 

u 

When  one  simulation  is  complete,  another  simulation  may  begin.  The 
random  number  generator  may  be  adjusted  to  reflect  the  natural  scatter  of 
a  given  group  of  test  specimens,  see  flow  diagram  in  Figure  2.9. 

* 

Simulation  Examples  -  Static  Loading. 

The  experimental  examples  shown  in  Figures  2. A  and  2.5  (quasi-static 
loading  case)  will  be  simulated  by  the  foregoing  discussed  stochastic 
search  routine.  The  laminates  considered  were  made  from  the  AS-3501-06 
graphite-epoxy  system.  The  nominal  single  ply  thickness  of  the  material 
is  0.005  in.  (0.125  mm).  The  unidirectional  ply  is  experimentally  char¬ 
acterized  as  one  of  monoclinic  symmetry  [4]  and  its  elastic  moduli  were 
reported  in  Ref.  [15]. 

In  the  calculation  Involving  equations  (14)  -  (20),  the  following 
laminate  properties  were  also  experimentally  characterized  [15]; 

&T  =  225°F  (125‘’C) 

t  =  0.0052  in.  (0.13  mm)  (21) 

2  2 

=  1.3  in-lb/in  (227.5  j/m  )  for  transverse  cracking. 

For  the  definition  of  the  "effective"  flaw  distributions  f(a)  and 

f(S)  in  the  90'’-layers  of  both  the  [0  /90  ]  and  [0„/90„]  laminates,  the 

z  1  s  Z  3  s 

it 

All  examples  in  this  chapter  are  from  the  test  data  reported  in 
Ref.  [15]. 


normal  distributions  as  expressed  in  (11)  and  (12)  were  adapted.  The 


parameters  V  and  o  were  chosen  as:  ^ 

v  =  0.0032";  a  =  0.0008" 

^  (0.08  mm)  ^  (0.02  mm) 

(22) 

V  =  0.0125";  a  =  0.0046"  * 

®  (0.313  nm)  ®  (0.115  mm) 

These  parameters  were  chosen  to  reflect  the  fact  that  the  largest 
flaw  in  the  set,  say  the  99  percentile  of  f (a)  ,  may  take  a  value:  I 

a  ~  p  +  3a  =  0.0056"  (23) 

(0.14  mm) 

and  that  for  =  0.0125",  it  implies  80  effective  flaws  for  each  inch  of  P 

the  90‘’-layer  along  the  length  of  the  specimen  (32  flaws  per  cm)  .  The 

choice  of  is  based  on  actual  test  data  from  sample  specimens. 

Having  defined  the  specific  form  of  the  effective  flaw  distribution,  * 

there  remain  three  other  inputs  of  information  in  order  to  start  the  simula¬ 
tion.  Namely,  the  calculated  energy  release  rate  function,  G^(a^,  a),  the 
calculated  strain  energy  retention  factor ,  R{S) ,  and  the  ultimate  laminate  ' 

tensile  strength ,  . 

For  a)  it  is  necessary  to  calculate  for  each  of  the  two 

laminates  their  respective  coefficient  functions  C^(a)  and  C^(a),  see 
Equation  (20).  Figures  2.10  and  2.11  show  the  calculated  C^(a)  and  C^(a), 
respectively. 

For  the  laminate  family  [0„/90„]  and  [0„/90„]  ,  a  single  R(S)  curve 

Z  Z  S  4  J  s 

can  be  generated  when  S  is  normalized  in  terms  of  the  90° -layer  thickness. 

This  is  shown  in  Figure  2.12.  The  R(S)  curve  can  be  fitted  by  a  simple 


exponential  equation  as, 


R(S/nt)  =  1  -  1.0056  exp  (-0.7776  S/nt)  .  (24) 

The  detail  of  the  R(S)  calculation  is  included  In  Appendix  A. 

Finally,  input  for  for  °u  ~  (958  mpa) 

and  for  ^u  ~  (821  mpa).  Both  were  suggested  by  ex¬ 

periment  . 

The  actual  simulation  was  carried  out  by  a  computer  routine  programmed 
on  a  VAX  11-750  computer. 

Figures  2.13  and  2.14  depict  the  transverse  cracking  simulation  results 
for  the  respectively,  under  quasi¬ 

static  tension.  In  each  case,  five  simulated  specimens  were  used.  The 
results  in  these  two  figures  are  to  be  compared  to  that  shown  respectively, 
in  Figues  2.4  and  2.5.  For  all  practical  purposes,  the  comparison  is 
extremely  close. 

More  simulated  transverse  cracking  results  have  been  generated  and 
compared  to  experiments  In  the  present  study.  These  are  found  in  Appendix 
B.  Generally,  the  simulation  model  describes  well  the  transverse  cracking 
process  under  static  tension. 


Figure  2.3  Schematics  of  Multiple  Transverse  Cracks  In  [0/90]^  Type 
Laminates. 
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Figure  2.5  Experimental  Transverse  Crack  Density  Vers\is  Applied  Laminate  Stress  3  . 
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Figure  2.9  Flow-Chart  For  Monte-Carlo  Simulation 
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2.4  Mechanics  of  Multiple  Transverse  Cracks  -  Fatigue  Loadln 


In  this  section,  the  energy  release  rate  concept  for  crack  growth 
will  be  extended  to  model  the  transverse  cracking  process  found  In  lami¬ 
nates  when  they  are  subjected  to  constant  amplitude  fatigue  loads. 

T^or  most  laminates,  fatigue  experiments  have  shown  that  the  develop¬ 
mental  process  of  transverse  cracks  Is  similar  to  that  found  under 
statically  applied  load.  The  effect  of  the  fatigue  loading  Is  essentially 
a  shift  of  the  damage  event  with  time  (or  fatigue  cycles).  Figures  2.15 
and  2.16  show  the  experimental  crack  density  versus  cycle  plots  for,  re¬ 
spectively,  under  different  levels  of 

fatigue  amplitudes.  Figure  2.17  Illustrates  the  similarity  between  the 
crack  development  processes  under  static  and  fatigue  loads.  The  simi¬ 
larity  between  them  suggests  the  existence  of  a  "damage"  function  which 
maps  the  damage  growth  curve  in  n  -  a  space  onto  n  -  N  space.  This 
mapping  function  relates  the  parameters  n,  o,  N  in  a  certain  form. 

Recall  in  the  previous  quasi-static  loading  case  that  the  effective 
flaws  are  regarded  as  the  aggregate  representation  of  the  inherent  material 
microflaws.  This  same  concept  may  be  equally  valid  in  the  case  of  fatigue 
loading.  Accordingly,  assume  here  that  the  sizes  {a^}  of  the  effective 
flaws  increase  continuously  with  time  during  fatigue  loading.  VJhen  the 
individual  size  reaches  a  characteristic  size  a*,  a  transverse  crack 
is  formed.  This  assumption  provides  the  necessary  criterion  for  a  flaw  to 
become  a  crack;  and  it  also  defined  the  time  under  load  when  the  crack  is 
formed  with  the  aid  of  a  flaw  growth  equation. 

The  previous  crack  slimilation  routine  can  thus  be  applied  to  fatigue 
Induced  transverse  cracking  such  as  shown  in  Figures  2.15  and  2.16, 


L 


provided  a  flaw  growth  law  is  formulated  for  the  "effective"  flaws  under 
cyclic  loads. 


Let  the  growth  of  the  effective  flaw  follow  the  conventional  power 
law  for  some  brittle  metals  [26] : 

j  G(a^,  a)  - 

H  =  a  [  -  ]P  (25) 

c 

where  a  and  p  remain  as  disposable  constants  and  is  the  amplitude  of 
the  applied  laminate  fatigue  load. 

Equation  (25)  contains  implicitly  5^,  a,  and  the  effects  of 
the  laminate's  layering  structure  through  the  quantity  G(5^,  a). 
Furthermore,  by  defining  G(a^,  a)  in  the  form  (19),  the  influence  of  exis¬ 
ting  transverse  cracks  on  a  flaw  is  also  Included.  The  material  resis¬ 
tance  against  flaw  growth  is  expressed  by  G^ .  The  parameters  ci  and  p  are 
expected  to  be  material  related  only. 

The  choice  of  the  parameters  a  and  p  in  (25)  will  clearly  involve 
some  degree  of  empiricisms.  However,  justification  for  the  choice  can 
come  only  from  a  careful  and  exhaustive  experimental  correlation.  It 
seems  that  a  verification  by  direct  physical  measurement  at  the  macroscopic 
dimensional  level  is  impractical  if  not  impossible. 

An  integration  of  (25)  gives: 


i  I  ,  -5,, 

a  ^c 


where  N  is  the  number  of  cycles  at  which  the  flaw  of  size  2a  becomes  a 
transverse  crack. 
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The  simulation  routine  outlined  in  Figure  2.9  is  slightly  modified 
to  model  the  transverse  cracking  process  under  fatigue  loads.  The  simu¬ 
lation  provides  the  desired  time-sequence  crack  formation  curve  under  the 
applied  cyclic  load. 

Simulation  Examples  -  Fatigue  Loading. 

In  the  simulation  of  fatigue  induced  cracks,  the  same  laminates 

[0_/90-]  and  [0.,/90.,]  will  be  considered.  The  criterion  for  the  forma- 
L  L  s  z  J  s 

tion  of  a  crack  from  a  flaw  of  size  2a  is  defined  when  a  =  a*,  a*  is 
set  equal  to  h  (half  the  thickness  of  the  90°-layer) .  The  law  which 
governs  the  growth  of  a  under  the  constant  amplitude  cyclic  load  takes 
the  form  of  (25),  with: 

a  =  0.004t;  p  =  20  (27) 

And,  of  course,  the  effective  flaw  distributions  f(a)  and  f(S)  in 

the  90'’-layer  remain  the  same  as  before  (see  Equations  (11)  and  (12)); 

their  parameters  have  the  same  value  as  expressed  in  Equation  (22).  In 

addition,  the  calculated  G^(oj,  a)  and  R(S)  functions  also  remain  the 

same  (see  Figures  2.10,  2.11,  and  2.12).  All  simulations  are  terminated 
*  6 

once  N  reaches  10  cycles. 

Figures  2.18  and  2.19  depict  the  transverse  cracking  simulation  re¬ 
sults  for  the  two  types  of  laminates  under  a  family  of  tensile  fatigue 
loads.  For  each  fatigue  load  amplitude,  three  simulated  specimens  were 
used.  These  results  are  to  be  compared  to  their  experimental  counterparts 
shown  in  Figures  2.15  and  2.16.  Again,  the  comparison  appears  quite  good. 
More  comparative  study  cases  between  the  simulations  and  experiments  for 
fatigue  induced  transverse  cracking  are  found  in  Appendix  C. 
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Figure  2.18  Simulated  Transverse  Crack  Density  Versus  Fatigue  Cycle  N.  Co 2/9^2^ 
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2.5.  A  Cumulative  Damage  Model  for  Multiple  Transverse  Cracks . 

From  the  previously  developed  model  for  multiple  transverse  crack 
growth  under  constant  amplitude  fatigue  load,  a  state  of  "constant 
damage"  can  be  quantitatively  defined.  For  instance,  the  trans¬ 
verse  crack  density  n  as  a  measure  of  the  "damage  state."  The  damage 
state  n  can  be  induced  either  by  the  constant  amplitude  fatigue  load 
in  N  cycles  or  by  the  static  load  o^.  It  is  assumed  that  a  unique 
0  -  N  relationship  could  be  determined  for  a  given  damage  state  n.  This 
relationship  can  be  found  by  means  of  the  previously  developed  fatigue 
crack  growth  model. 

Figure  2.20  Illustrates  schematically  a  family  of  constant  damage 
lines  in  the  o  -  N  space.  These  curves  are  calculated  numerically  from 
the  fatigue  crack  growth  simulation  routine  by  setting  n  at  some  given 
constant  value.  Note  that  these  constant  damage  curves  can  also  be  gen¬ 
erated  from  the  experimental  crack  density  growth  curves. 

From  the  constant  damage  curves  in  the  a  -  N  plane,  a  cumulative 
damage  model  for  the  crack  growth  process  under  spectrum  load  is  construc¬ 
ted  as  follows. 

Consider  the  spectrum  load  history  such  as  shown  in  Figure  2.21(a). 
The  laminate  is  subjected  first  to  a  constant  amplitude  fatigue  load 
for  a  period  of  N^^  cycles.  The  damage  at  the  end  of  cycles  is  found 

at  point  A  in  the  a  -  N  domain.  Figure  2.21(b).  The  constant  damage 
curve  passing  point  A  is  denoted  by  Then,  the  amplitude  of  the 

fatigue  load  is  increased  to  02  for  a  period  of  (N2  -  N^)  cycles,  see 
Figure  2.21(a).  Since  the  damage  in  the  laminate  at  the  onset  of  02  is 
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still  N^,  the  location  of  this  damage  state  is  at  the  intersect  of  02 
and  the  curve  point  B.  For  a  period  of  (N^  -  N^)  fatigue  cycles, 
the  damage  state  in  the  laminate  advances  to  C,  which  is  located  on  the 
constant  damage  line  N2.  At  the  end  of  N2  cycles,  the  fatigue  amplitude 
is  lowered  to  again  for  a  period  of  (N^  -  N2)  cycles.  During  this 
period,  the  damage  state  advances  from  point  D  on  curve  N2  to  point  E 
on  curve  . 

At  each  level  of  the  fatigue  load,  the  initial  damage  state  is  de¬ 
termined  if  the  number  of  cycles  is  specified.  Actually,  the  damage 
state  at  the  end  of  the  spectrum  loading  is  found  from  integrating  the 
damage  growth  equation: 
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(28) 


This  method  can  easily  be  extended,  on  a  cycle  by  cycle  basis,  to 
random  spectrum  loadings. 


Simulation  Examples  -  Spectrum  Loading. 

Let  us  consider  again  the  laminates  [0„/90„]  and  [0„/90_]  .  The 

2  2s  2  i  s 

transverse  crack  growth  processes  in  these  laminates  under  constant  ampli¬ 
tude  fatigue  loading  have  been  simulated  by  the  proposed  stochastic  pro¬ 
cedure.  Families  of  crack  density  (n)  versus  fatigue  cycles  (N)  were 
generated  and  shown  in  Figures  2.18  and  2.19.  From  these  crack  density 
growth  curves  in  the  n  -  N  space,  a  family  of  constant  damage  lines  may  be 
obtained  by  a  graphical  method.  This  is  done  by  drawing  a  constant  n  line 
in,  say.  Figure  2.18.  The  intersecting  points  with  the  n(N)  family  curves 
yield  the  desired  constant  n  curve  in  the  a  -  N  space.  Figure  2.22  shows 
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the  constant  damage  family  lines  for  Figure  2.23  shows 

the  same  for  [0„/90„]  . 

i  J  s 

The  constant  damage  lines  will  be  used  to  construct  the  desired 
cumulative  damage  model  which  predicts  the  transverse  cracking  state  in 
the  laminates  under  a  given  variable  amplitude  (spectrum)  fatigue  load. 
The  procedure  has  already  been  explained  earlier  in  this  section. 

To  illustrate  the  accuracy  of  the  cumulative  model,  a  number  of 
variable  loading  cases  have  been  considered.  Four  typical  test  cases 
between  the  predicted  results  and  the  exerpimental  results  are  displayed 
in  the  following  table.  In  each  case,  the  prediction  is  based  on  the 
graphical  determination  such  as  illustrated  in  Figure  2.21.  The  graphs 
used  in  the  four  test  cases  are  shown  in  Figures  2.23  and  2.24.  More 
theoretical/experimental  comparison  may  be  found  in  Appendix  D  of  this 
report . 


Generally,  the  comparison  is  very  close. 
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Figure  2.21  Schematics  ior  the  Concept  of  a  Cumulative  Dawa^e  Model 
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2.6.  Mechanics  of  Fre e  Edge  De lamination  -  Static  Loading . 

Free  edge  induced  delamination  crack  growth  has  been  frequently 
observed  as  the  most  damaging  sub-laminate  failure  mode.  The  crack  de¬ 
velopment  process  is  generally  idealized  as  a  plans  crack  which  is  formed 
along  the  laminate  free  edge,  propagating  inward  along  an  interface  of 
two  adjacent  layers.  But,  post-failure  examination  of  free  edge  delami¬ 
nation  of  ten  showed  a  complicated  crack  growth  path;  the  path  is  usually 
not  contained  in  one  given  layer  interface.  For  example,  consider  a 
laminate  in  the  form,  say,  [+25/90]^  under  uniaxial  tension.  Because  the 
laminate  contains  a  90'’-layer,  transverse  cracking  in  this  layer  could  be 
Induced  by  the  applied  tensile  load.  However,  the  thickness  of  the  90°- 
layer  is  not  large;  a  relatively  high  load  is  required  to  produce  the 
transverse  crack,  according  to  the  energy  theory  presented  earlier.  In 
this  case,  then,  a  free  edge  delamination  may  be  induced  before  transverse 
cracking.  Once  edge  delamination  is  started,  multiple  transverse  cracks 
in  the  90°-layer  could  still  be  formed  inside  the  delamination  area.  Fre¬ 
quently,  the  delamination  growth  process  is  complicated  by  the  presence  of 
transverse  cracks,  e.g.,  the  delamination  plane  can  be  shifted  to  the 
25/90  Interface  because  of  the  transverse  cracks  in  the  90°-layer. 

As  a  counterexample,  consider  the  laminate  of  [+25/90-]  ,  also  under 
uniaxial  tension.  In  t.iis  case,  transverse  cracking  will  he  the  first  to 
occur  because  of  the  larger  90'’-layer  thickness.  The  presence  of  the 
transverse  cracks,  in  turn,  causes  delamination  In  the  25/90  interface. 

It  seems  that  the  actual  thickness  of  the  90°-layer  in  the  laminate  can 
alter  the  competing  damage  modes  and  their  occurrence  sequences.  In  this 
section,  we  shall  confine  our  discussion  to  the  mechanisms  of  delamination 
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alone,  not  complicated  by  other  interacting  failure  modes.  The  latter 
subject,  though  extremely  Important,  will  be  studied  separately. 

Analytical  studies  of  the  free  edge  delamination  problem  have  origi¬ 
nated  from  the  much  celebrated  work  of  Pipes  and  Pagano  [28],  who  formu¬ 
lated  and  computed  the  boundary-layer  interlaminar  stress  solutions  for 
a  long  symmetric  laminate  under  tension.  The  formulation  is  based  on 
the  macroscopic  ply-elasticity  theory,  which  regards  the  material  layers 
as  individually  homogeneous  media;  material  and  geometrical  discontinu¬ 
ities  exist  only  across  the  layer  interfaces. 

Generally,  the  free  edge  stress  field  is  three-dimensional  and  is 
singular  at  the  Intersection  of  the  free-edge  with  the  layer  interface  [29]. 
Clearly,  interface  delamination  is  caused  by  the  highly  concentrated 
edge  stresses,  especially  the  interlaminar  normal  stress  0^. 

In  a  series  of  tensile  strength  tests,  Bjeletich,  et  al.  [30]  examined 
the  failure  modes  of  six  families  of  quasi- isotropic  laminates  by  alter¬ 
nating  the  stacking  sequence  of  the  0°,  90°  and  +45°  layers.  Edge  stress 

analysis  indicated  a  compressive  a  along  the  free  edge  of  the  [0/90/+45] 

z  s 

laminate,  while  a  tensile  a  for  the  [+45/0/90]  laminate.  The  latter 

z  —  s 

developed  premature  delamination  and  the  growth  of  it  had  led  to  a  much 
lowered  tensile  strength.  Though  the  knowledge  of  the  free-edgt  stresses 
can  provide  an  explanation  why  delamination  occurs,  a  quantitative  pre¬ 
diction  for  its  occurrence  requires  an  independently  defined  criterion. 

In  a  study  on  delamination  for  similar  graphite-epoxy  laminates, 
Rodini  et  al.,  [31]  determined  experimentally  the  critical  tensile  loads 
at  the  onset  of  free  edge  delamination  in  a  ~  ^ 

family.  They  found  that  the  critical  laminate  tensile  stress,  o^,  varied 


greatly  with  the  value  of  n.  Specifically,  the  critical  decreases 

at  the  rate  of  about  /n,  even  though  an  edge  stress  analysis  yields 

Identical  a  for  the  same  a  for  all  values  of  n. 
z  X 

As  we  now  know  more  about  the  failure  sequence  In  laminates  such 
as  [+45/0/90]^,  the  results  of  Rodini  et  al.,  [31]  evidently  also  in¬ 
volved  transverse  cracking  In  the  90“-layers.  It  is  possible  that  all 
three  laminates  (n  =  1,  2,  3),  or  at  least  those  of  n  =  2  and  3,  had 
transverse  cracks  before  delamlnatlon.  Hence,  in  order  to  analytically 
describe  the  delamlnatlon  process  generally,  one  must  take  Into  account 
the  interacting  effects  of  the  transverse  cracks. 

The  Energy  Criterion 

Wang  and  Crossman  [13]  assumed  that  material  flaws  also  exist 
randomly  on  any  one  of  the  Interfaces  between  the  material  layers. 

Those  flaws  which  are  located  within  the  free-edge  stress  zone  from  an 
"effective"  flaw  having  a  size  a^  at  the  instant  of  onset  delami¬ 
nation,  see  inset  in  Figure  2.24.  Note  that  a^  is  still  an  unknown 
quantity.  And,  it  can  also  be  random  in  nature. 

Following  the  finite  element  crack-closure  procedure  by  allowing 
virtual  extension  of  a  along  the  layer  interface,  energy  release  rate 
curves  for  externally  applied  load  are  generated  (for  detail,  see  [13]). 

Figure  2.24  shows  a  typical  G  curve.  The  energy  release  rate  G 
increases  sharply  with  a,  but  reaches  an  asymptote  at  a  =  a^.  Generally, 
the  value  of  a  is  about  one-half  the  layer-thickness  which  contains  the 
delamination  [13].  Thus,  the  layer  thickness  may  in  some  cases  affect 


the  value  of  a  ,  and  also  the  value  of  G. 


Expressing  the  computed  energy  release  rate  in  the  form 


G(e  ,  a)  =  [  /C  (a)  •  e  +  /C^(a)  -AT]  •  t 

o  e  o  T 


(29) 


and  applying  the  Griffith  criterion  (14)  for  the  onset  of  the  delamination 
crack,  one  obtains  the  critical  far-field  strain; 


/G  /t  -  /C_(a)  •  AT 
c _ i _ 

/C  (a) 


a  a  . 
o 


(30) 


The  problem  comes  dovm  to  two  important  questions;  namely,  what  value 

is  a  and  in  which  interface  is  it  located? 
o 

To  answer  the  first  question,  recall  our  earlier  discussions  about 
the  effective  flaw  distribution  f(a)  in  the  90®-layers.  The  worst  flaw 
may  have  a  size  equal  to  one  ply  thickness  for  commercial  graphite-epoxy 
systems.  It  is  believed  that  the  worst  flaw  which  is  responsible  for 
delamination  is  at  least  of  this  magnitude  if  not  larger,  because  of 
possible  additional  cutting  flaws  along  the  free  edge.  In  any  event,  if 
a^  in  Figue  2.24  is  less  than  2  times  the  ply  thickness,  one  simply 
uses  a  =  a^  and  predicts  from  Equation  (30)  the  minimum  possible  load  for 
the  onset  of  delamination.  Thus 
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The  answer  to  the  second  question,  however,  requires  additional  analy¬ 
sis.  Consider  as  an  example  a  ^^2^s  uniform  compression. 

Interlaminar  tensile  is  developed  in  this  laminate  and  edge  delamina¬ 
tion  is  Induced  before  global  buckling,  see  Figure  2.25.  A  three-dimen¬ 
sional  through-thickness  display  of  is  shown  in  Figure  2.26.  The 

largest  O  is  located  on  the  45/-45  interface  (in  fact,  this  a  is 
z  ’  z 


singular  at  the  free  edge) .  On  the  other  hand,  a  along  the  mld- 

z 

plane  (z  =  0)  is  finite  at  the  free  edge.  In  addition,  there  is  also 
a  singular  interlaminar  shear  stress  along  the  45/-45  interface. 

Thus,  the  analysis  suggests  the  0/45  interface  as  the  most  probable 
delamination  site. 

However,  according  to  the  energy  release  rate  argument  as  well  as 
the  "effective"  flaw  hypothesis,  one  must  rely  upon  the  calculated  energy 
release  rate  curve  in  order  to  define  which  one  of  the  interfaces  is 
likely  to  delaminate. 

The  finite  element  crack-closure  results  of  the  G(a)  curves  for  this 
example  problem  is  shown  in  Figure  2.26.  Indeed,  the  interface  that 
yields  the  highest  possible  energy  release  rate  is  the  mid-plane  (z  =  0) , 
not  the  45/-45  interface  (though  the  reverse  is  true  when  a  is  smaller 
than  one  ply  thickness).  Besides,  a  delamination  in  the  mid-plane  is  in 
mode-I,  while  it  is  mixed-mode  (I,  III)  in  the  45/-45  interface.  Mixed- 
mode  is  usually  higher  than  Thus,  the  mid-plane  is  the  most  prob¬ 

able  delamination  site,  not  the  0/45  interface. 


Numerical  Examples  -  Static  Loading. 

In  order  to  assess  the  accuracy  of  the  energy  model  presented  above, 
experiments  have  been  conducted  using  AS-3501-06  graphite-epoxy  [+25/90^]^, 
n  =  1,  2,  3  (see  Reference  [15],  for  experimental  details).  Onset  and 
growth  of  free  edge  delamination  were  recorded  in  all  cases.  Figure  2.27 
displays  the  delamination  versus  the  applied  load  plots.  From  these  plots, 
the  average  onset  load  for  each  laminate  case  can  be  extrapolated.  The 
growth  of  delamination  is  essentially  stable,  but  the  growth  pace  is 
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rather  rapid.  In  all  cases,  delamlnatlon  growth  led  to  final  failure. 

Delaminatlon  in  [+25/90-]  and  [+25/90_]  started  after  the  occurrence 

2  s  j  s 

of  transverse  cracks  in  the  90°-layer.  Hence,  we  shall  defer  the  discuss¬ 
ion  on  these  two  laminates  to  later  in  this  section.  Here,  we  shall  con¬ 
sider  onlv  [+25/90]  . 

s 

A  free  edge  stress  analysis  and  an  energy  release  rate  calculation 
Indicate  that  for  the  [+25/90]^  laminate,  mid-plane  delamination  is  most 
probable.  The  energy-release  rate  curves  and  as  a  function  of 
delamination  size  a  in  the  mid-plane  are  displayed  in  Figures  2.28  and 
2.29,  respectively. 

2 

Using  the  energy  curves  and  the  value  of  G  =  1.0  in-lb/in  for  mode-I 

c 

delaminatlon  growth  (see  discussion  on  in  Section  2.2),  the  minimum 
load  at  the  onset  of  edge  delamination  is  calculated  via  Equation  (31) . 

The  following  gives  a  comparison  between  the  predicted  and  the  experi¬ 
mental  results  for  the  delaminatlon  onset  load: 

Experimental,  Ksi  Predicted,  Ksi 

48  -  50  (3  samples)  48  (Average) 

The  prediction  is  based  on  a  mid-plane  delamination  while  the  experi¬ 
mental  observation  of  the  delaminated  specimens  Indicated  a  delaminatlon 
path  which  zig-zags  between  the  90/90  and  the  25/90  interfaces.  Even  with 
these  practical  differences  between  the  theory  and  the  experiment,  the 
energy  model  has  yielded  exceedingly  satisfactory  results.  More  results 
from  this  study  are  found  in  Appendix  E. 
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Figure  2.26  Interlaminar  Normal  Stress  Distribution  and  Energy  Release 
Rate  Coefficients  for  Compression. 
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Figure  2.27  Delamlnatlon  Growth  (a/w)  versus  Applied  Load  for  (+25/90^ 


Figure  2.29  Enegry  Release  Rate  Coefficients  For  Delamination 
Under  Thermal  Load. 


It  is  well  known  that  cyclic  fatigue  loads  also  induce  sublaminate 
matrix  cracks.  Including  many  forms  of  delamination  growth.  In  fact,  such 
cracks  can  be  propagated  at  fatigue  stress  levels  much  below  the  threshold 
stress  determined  statically.  Laboratory  tests  tend  to  indicate  that  fatigue 
load  induced  damage  process  is  generally  similar  to  that  induced  by  static 
load;  and  the  effects  of  cyclic  loading /unloading  are  manifested  by  a  shift- 
of-time. 


Figures  2.30,  2.31,  and  2.32  show  the  experimental  delamination  (in  terms 
of  percent  specimen  width)  versus  the  number  of  cycle  plots  for  the  [+25/90]^ 
laminate  under  the  fatigue  stress  of  29,  33,  and  40  ksi,  respectively. 
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The  crack  growth  rate  law  which  is  used  widely  for  fine-grain  metallic 
material  is 


=  5  [G(a^,  AT,  a)/Gj‘ 


where  is  the  fatigue  stress  amplitude,  N  is  fatigue  cycle  and  a  and  p  are 

two  undetermined  empirical  constants.  ? 

Since  in  Equation  (32)  the  quantity  G(a^,  AT,  a)  includes  already  the 
laminate  stacking  sequence,  lamina  thicknesses,  laminate  shape  features,  the 
crack  geometry,  and  the  nature  of  loading  and  since  the  quantity  G/G^  repre¬ 
sents  the  relative  magnitude  of  the  crack  driving  force,  then  a  and  p  are 
expected  to  be  pure  material  constants.  And,  in  particular,  the  exponent  p 
is  the  most  Important  measure  of  the  crack  growth  rate. 

Recall  that  the  same  crack  growth  rate  law  (32)  was  used  in  the  stochas¬ 
tic  simulation  of  multiple  transverse  cracks  in  [0/90]  type  of  laminates 

“  I 
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stant  a  =  2  X  10  inch  and  p  =  20.0  were  used  for  mode-I  growth. 

At  this  point,  it  is  perhaps  appropriate  to  review  some  other  related 
works.  Wilkins  [19]  conducted  fatigue  growth  experiments  on  mode-I  and 
mixed-mode  delamination.  Th>_  test  results  suggest  that  a  power  law  such  as 
Equation  (32)  would  govern  the  delamination  growth.  For,  if  taking  logarithm 
on  both  sides  of  (32),  a  linear  relationship  between  da/dN  and  G(o^,  a)  is 
obtained;  the  constant  ratio  between  them  -lelds  the  exponent  p.  Wilkins' 
results  confirmed  the  adequacy  of  the  crack  growth  law  in  the  form  of  (32), 
for  both  mode-I  and  mixed-mode  growth. 

But,  Wilkins  also  pointed  out  that  the  exponent  p  for  mode-I  growth  is 
much  larger  than  that  for  mixed-mode  growth.  For  the  material  used  in  his 
study  (AS-3501-06) ,  he  found  p  ~  20  for  mode-I  and  p  -  8  for  mixed-mode  de- 
lamination  growth.  Since  the  value  of  G/G^  in  (32)  is  normally  smaller  than 
unity,  a  larger  p  (positive)  Implies  slower  growth  rate.  It  has  been 
suggested  that  mode-I  delamination  is,  generally,  a  static  issue;  and  mixed¬ 
mode  delamination  is  a  fatigue  issue,  so  far  as  epoxy-based  composite 
laminates  are  concerned. 

To  fit  the  experimental  data  shown  in  Figures  2.30,  2.31  and  2.32,  the 

corresponding  mixed-mode  energy  release  rate  curves  are  used,  see  Figures 

2 

2.28  and  2.29.  And  for  mixed-mode  cracking  G^  =  1.3  ar.-lb/in  .  The  value 
of  a  and  p  are  set  at 

a  =  0.04  inch;  p  =  8.3  (33) 

Thus,  the  predicted  delamination  growth  curves  as  a  function  of  fatigue 
amplitude  and  fatigue  cycles  are  shown  in  solid  lines  in  Figures  2.30,  2.31 


Figure  2.30  I'ercunt  UelamtnaC Ion  vs.  Fatigue  Cycles  for 
[+25/90 ]g  Laminates,  0“29  KSl. 


lamination 


Figure  2.31  Percent  Delai.iination  vs.  Fatigue  Cycles  tor 
[+25/90Jg  Laminates,  o=33  KSI. 


2.12  Dc I  .iiiii ii.k  inn  vs.  F.ir  i  (!yt  l 

(+25/90 1 j.  Laminates,  o-AO  KSI. 


and  2.32.  More  comparative  results  on  this  particular  laminate  under 
other  load  amplitudes  are  Included  In  Appendix  F. 

2.8.  A  Cumulative  Damage  Model  for  Free  Edge  Delamlnatlon 

As  shown  In  the  previous  section,  the  growth  of  edge  delamlnatlon  can 
be  satisfactorily  simulated  by  a  suitable  crack  growth  law.  Specifically, 
Equation  (32)  Is  written  In  a  form  which  contains  the  loading  factor, 
lamination  geometry,  the  basic  material  moduli,  etc.,  through  the  computed 
quantity  GCc^,  a).  The  previously  developed  constant  damage  concept  de¬ 
scribes  the  delamlnatlon  growth  accumulation  process  for  the  same  laminate 
under  variable  amplitude  fatigue  loads.  The  procedure  Is  similar  to  that 
used  for  the  transverse  cracking  problems  discussed  In  Section  2.5.  Upon 
Integrating  Equation  (32),  a  family  of  constant  damage  lines  can  be  gener¬ 
ated  numerically  or  graphically  In  the  o  -  N  space.  Take,  as  an  example  the 
[+25/90]^  laminate's  results  shown  In  Figures  2.30,  2,31,  and  2.32;  a  set  of 
constant  damage  lines  can  be  obtained  using  these  curves.  The  latter  Is 
shown  In  Figure  2.33  where  the  state  of  damage  Is  expressed  In  terms  of  the 
normalized  size  of  delamlnatlon,  a/w  (w  Is  the  half-width  of  the  specimen) . 
The  constant  damage  curves  can  predict  the  amount  of  cumulative  damage 

(delamlnatlon  size)  In  the  laminate  under  spectrum  loads. 

I 

As  examples,  consider  first  a  case  where  the  [+25/90]  laminate  Is 

s 

under  Oe  “  20  ksl  (138  rapa)  for  10^  cycles  and  Is  then  loaded  by  = 

I  30  ksl  (207  mpa)  for  3  x  10^  cycles.  To  predict  the  cumulative  damage,  we 

used  the  curves  In  Figure  2.33.  Under  ’*  20  ksl  for  10^  cycles,  there 
will  be  no  delamlnatlon  according  to  the  constant  damage  lines.  When  the 

! 
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load  is  increased  to  30  ksi,  a/w  =  25%  is  predicted.  The  experi¬ 
mental  results  for  this  case  show  a  21%  delamination  at  the  end  of 
loading,  see  Figure  2.33. 

Another  case  concerns  the  same  laminate  loaded  under  =  20  ksi 
for  3  X  10^  cycles,  followed  by  cy£  *  ^0  ksi  (276  mpa)  for  500  cycles, 
followed  by  =  30  ksi  for  2  x  10^  cycles.  At  the  end  of  each  loading 
step,  the  predicted  damage  state  is,  a/w  =  0,  3%  and  25%,  respectively, 
while  the  experiment  indicated  a/w  =  0,  4%  and  33%,  see  Figure  2.34. 

The  experimental  results  for  Che  above  tvo  cases  have  been  reported 
in  Ref.  [15].  A  more  detailed  comparison  between  the  prediction  and  the 
experiment  for  these  two  cases  is  given  below; 


Fatigue  Amplitude 

Ksi  (mpa) 

Cycles  @  Load 

Delamlnatlon,  %w 

Experiment  Prediction 

20 

(173) 

10^ 

0 

0 

30 

(259) 

10^ 

10 

12 

3  X  10^ 

21 

25 

20 

(173) 

3  X  10^ 

0 

0 

40 

(346) 

10^ 

1 

1 

5  X  10^ 

4 

3 

30 

(259) 

10^ 

25 

20 

2  X  10^ 

33 

25 

The  above  shows  that  the  prediction 

damage  lines  also  appears  satisfactory. 

based  on 

Many  more 

the  constant 

comparative  cases  have 

been  conducted  in  this  research.  All  the  results  are  found  in  Appendix  G. 
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2.9.  Edge  Delaminatlon  and  Transverse  Cracking  Interac t Ion s 

In  the  previous  examples,  the  laminates  considered  are  such  that 
edge  delamination  is  the  first  and  only  sub-laminate  failure  mode;  the 
onset  load  for  delamination  could  be  determined  accurately  using  the 
energy  model.  In  most  other  cases,  however,  edge  delamination  and  trans¬ 
verse  cracking  are  competing  failure  modes  which  interact  with  each 
other  through  a  complicated  mechanism. 

Recall  the  laminate  [+25/90]  which  was  discussed  earlier.  The 

—  s 

thickness  of  the  90°-layer  is  made  small  so  tliat  transverse  cracking  is 
not  possible  before  edge  delaminatlon.  Nonetheless,  the  stress  field  in 
the  90°-layer,  which  contains  the  delamination,  is  actually  very  complex, 
so  much  so  that  high  density  transverse  cracks  could  be  observed  during 
the  entire  delamination  process. 

On  the  other  hand,  in  the  tests  of  the  [+25/90^]^  laminate,  the  lami¬ 
nate  suffered  first  transverse  cracks  in  the  90°-layer;  subsequently,  edge 
delamination  in  the  25/90  interface  was  developed  at  a  slightly  higher 
load.  And,  the  delaminatlon,  in  turn,  caused  many  moie  transverse  cracks. 
Similar  failure  sequence  was  also  observed  in  the  [+25/90^]^  laminates. 
These  examples  clearly  demonstrate  the  profound  interacting  effects 
among  competing  failure  modes,  which  exist  in  most  practical  laminates. 

In  an  attempt  to  simulate  the  problem  approximately.  Law  [32]  assumed 

that  when  transverse  cracks  occurred  in  the  90'’-layer,  the  tensile  modulus, 

E_,  and  the  Poisson  ratio,  V  of  the  90“-layer  will  both  reduce  to  nearly 
T  i  L 

zero.  With  this  reduced  moduli  in  the  90'’-layer,  the  overall  laminate 
becomes  energetically  delamination  prone  in  the  25/90  interface,  resulting 


P 


P 


9 


P 
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In  a  mixed-mode  cracking  (I,  II).  Then,  by  simulating  a  mixed-mode 

delamination  in  the  25/90  interface  (for  details  see  [32]),  the  total  - 

energy  release  rate  G  =  G^.  +  G^^  can  be  calculated  using  the  same  two- 

dimensional  crack-closure  routine.  For  instance,  the  computed  and 

C  curves  for  [+25/90-]  and  [+25/90_]  are  shown  in  Figures  2.28  and 
1  —  zs—  Js  I 

2 

2.29.  Using  these  curves  and  setting  G^  =  1.0  in-lb/in  for  mixed-mode 
delamination,  the  minimum  onset  load  for  edge  delamination  can  be  pre¬ 
dicted.  A  comparison  between  the  prediction  and  the  experiment  is  given  - 

as  follows: 

Onset  Laminate  Stress.  Ksl  (mpa) 

Laminate  Experiment  Prediction 

[+25/9O2]  40  -  43  41  (283) 

(+25/90^1  32  -  33  38  (262) 

This  shows  that  the  prediction  for  [+25/90^]^  was  not  satisfactory.  ' 

The  short-fall  in  this  case  may  stem  from  the  fact  that  the  interaction 
effect  could  not  be  adequately  addressed  by  the  assumption  made  by  Law 
[33].  Further  refinement  of  the  simulation  model,  which  includes  failure 
modes  Interaction,  remains  to  be  developed. 

In  fact,  such  a  model  development  has  been  planned  in  the  Phase-III 
study  of  this  research  program.  At  this  reporting  time,  however,  no 
simulation  results  have  been  available  for  a  comparison  with  the  experi¬ 
mental  results  obtained  for  [+25/90„]  and  [+25/90-]  laminates. 

Z  s  3  s 
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III.  EXPERIMENTAL  PROGRAM  AND  RESULTS 


3.1.  Scope  of  Experiment  of  Phase-I. 

In  Phase-I  of  the  present  research  project,  the  experimental  program 
consisted  of  four  major  tests: 

(a)  material  and  laminate  characterization  tests; 

(b)  static  load-damage  growth  tests; 

(c)  constant  amplitude  fatigue-damage  growth  tests;  and 

(d)  cumulative  damage  growth  under  variable  amplitude  fatige 
load  tests. 

The  scope  of  test  (a)  was  to  characterize  the  basic  material  and 
laminate  properties  in  order  to  qualify  the  material  system  used  as 
well  as  to  obtain  the  basic  material  constants  needed  for  stress  and 
fracture  analyses. 

The  scope  of  test  (b)  was  to  identify  the  various  failure  modes  in 

laminates  under  static  loading  conditions.  The  tests  were  designed  to 

yield  Information  about  the  onset,  growth,  and  final  failure  mechanisms 

of  transverse  cracking  and  free  edge  delamlnatlon.  Essentially,  the 

laminates  [0„/90„]  ,  [0-90/_]  ,  and  [+25/90  ]  ,  n  =  1,  2,  3  were  tested. 
^  /  s  /  3s  —  ns 

These  laminates  were  also  tested  in  (c)  and  (d) ;  data  from  these 
tests  were  used  to  correlate  the  fatigue  crack  growth  equations  and  the 
cumulative  damage  model. 

All  test  results  obtained  in  Phase-I  were  collectively  reported  in 
the  Phase-I  final  report  [15].  Some  of  these  results  were  also  used  in 
Chapter  II  of  this  report  in  the  presentation  of  the  various  failure 


A  complete  presentation  of  the  test  data  along  with  the  analytical 
results  can  also  be  found  in  Appendices  A  through  H. 

The  details  of  the  test  methods  and  other  related  information  are 
contained  in  Ref.  [15]. 


.  Scope  of  Experiment  of  Phase-II 

In  Phase-II  of  this  research,  new  specimens  had  to  be  made.  Though 
using  the  same  AS-3501-06  material  system  as  in  Phase-I,  the  new  batch 
of  specimens,  however,  were  fabricated  at  the  University  of  Delaware 
(while  specimens  used  In  Phase-I  were  fabricated  at  NADC) .  For  this 
reason,  key  material  and  laminate  characterization  tests  had  to  be  per¬ 
formed  again  in  Phase-II. 

The  teat  details  were  malntaned  the  same  as  those  reported  in  Ref . 
[15].  A  complete  account  of  the  test  program  is  given  below: 

Laminate  Purpose  of  Test  Load  Type  Specimens  Used 

[Og]  Basic  Data  Static  to  faxi  18 

[90^2]  Basic  Data  Static  to  Fail  2 

(0.,/90„]  Basic  Data  Static  to  Fail  2 

l  Z  s 

T.C.  Growth  Step  Load  8 

T.C.  Growth  Fatigue  Load  2o 

(5  Fatigue  Loads) 

T.C.  Growth  Variable  Fatigue  ^ 

(3  cases) 


102/903!^ 

Basic  Data 

T.C.  Growth 

T.C.  Growth 

T.C.  Growth 

Static  to  Fail 

Step  Load 

Fatigue  Load 
(5  Fatigue  Loads) 

Variable  Fatigue 
(3  cases) 

2 

6 

17 

6 

[+25/90]^ 

Basic  Data 

Static  to  Fail 

2 

Edge  Delamination 

Static  to  Fail 

5 

Edge  Delamlnatlon 

Fatigue  Load 
(6  Fatigue  Loads) 

20 

Edge  Delamlnatlon 

Variable  Fatigue 
(5  cases) 

11 
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Laminate 

Purpose  of  Test 

Load  Type 

Specimens  Used 

[+25/90 J 
—  2  s 

Basic  Data 

Static  to  Fail 

2 

Interaction 

Step  Load 

5 

Interaction 

Fatigue  Load 

C. 

(2  Fatigue  Loads) 

[+25/90.] 

*“  J  s 

Basic  Data 

Static  to  Fail 

2 

Interaction 

Step  Load 

5 

Interaction 

Fatigue  Load 

(2  Fatigue  Loads) 

TOTAL  SPECIMENS  151 

All  data  obtained  in  Phase-II  are  presented  in  Appendices  A  through  H. 

Except  for  the  [+25/90.]  and  [+25/90_]  laminates,  all  the  data  has  been 
—  ^  s  —  J  s 

used  in  correlation  with  the  analytical  model.  Results  from  the  model  are 
listed  in  these  Appendices,  along  with  their  experimental  counterparts. 
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]V.  CURRENT  STATUS  OF  RESEARCH  AND  METHODOLOGY  SUMMARY 
4.1.  Current  Status 

The  principal  objective  of  this  research  has  been  illustrated  and 
discussed  in  Chapter  I.  The  details  of  the  model  development  are  presented 
in  Chapter  II.  At  this  reporting  time,  the  status  of  the  research  project 
may  be  summarized  as  follows: 

(a)  A  concept  of  "effective  flaw  distribution"  has  been  introduced. 
This  concept  postulates  that  natural  flaws  in  the  basic  material 
ply  (unidirectional)  and/or  between  ply  interface  (interlaminar 
flaws)  could  be  effectively  represented  by  a  suitable  distribu¬ 
tion  function.  This  concept  Is  then  used  in  conjunction  with 
energy  method  of  the  classical  fracture  mechanics.  Together 
they  provide  a  methodology  for  mu'^tiple  failure  sequences  in 
the  laminate. 

(b)  A  simple  crack  growth  law,  based  on  the  energy  argument,  has 
been  introduced  and  describes  the  individual  effective  flaw 
growth  behavior  under  constant  amplitude  fatigue  loading.  This 
crack  growth  law  implicitly  or  explicitly  contains  the  material 
factor  (fracture  toughness),  the  load  factor  (nature  of  load, 
amplitude,  etc.)  and  the  geometrical  factor  (ply  thickness, 
stacking  sequence,  free  edges,  etc.)  of  the  considered  laminate. 

(c)  The  concept  of  constant  damage  state  has  been  introduced.  It 
describes,  in  quantitative  terms,  the  state  of  damage  in  load- 
cycle  (a  -  N)  space.  This  concept  is  used  to  form  a  simple 
cumulative  damage  model  when  the  laminate  is  subjected  to 
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variable  amplitude  fatigue  load. 

(d)  All  the  above  have  been  applied  to  the  failure  modes  of  trans¬ 
verse  cracking  and  free  edge  delamlnatlon.  The  modelling  of 
these  two  failure  modes  has  been  successfully  correlated  with 
the  experimental  test  results. 

(e)  The  presently  developed  methodology  has  been  applied  to  simple 
laminates,  where  only  one  single  failure  mode  Is  Induced  by  the 
load.  It  remains  to  demonstrate  that  the  modelling  method  can 
be  applied  equally  well  to  laminates  having  complicated  lamina¬ 
tion  geometry.  Complication  may  arise,  however,  when  two  or 
more  failure  modes  occur  simultaneously  and/or  sequentially. 

In  these  cases,  modes'  Interaction  will  be  encountered.  We  feel 
that  the  methodology  should  serve  well  for  these  cases.  The 
necessary  numerical  tool  must  be  developed  to  handle  the  expected 
three-dimensional  stress  fields  when  two  or  more  Interacting 
singularities  are  present. 

(f)  For  laminates  which  suffer  extensive  delamlnatlon  before  final 
failure.  It  seems  practical  to  use  the  "delamlnatlon  state"  as 
a  criterion  for  final  failure.  However,  this  criterion  is  not 
expected  to  be  a  general  rule,  as  in  many  laminates,  delamlnatlon 
may  not  govern.  This  is  especially  true  for  laminates  which  con¬ 
tain  a  high  percentage  of  O'-plies.  Clearly,  it  is  still  a  long 
way  before  we  understand  the  complexity  of  the  damage  processes 
when  the  laminate  Is  approaching  final  failure.  Until  then,  the 
development  of  a  general  laminate  strength  criterion  (Including 
the  S  -  N  curve  in  the  o  -  N  space)  can  only  be  ad-hoc  in  nature. 
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4.2.  Methodology  Sunmary 

This  sunmary  Is  Intended  to  illustrate  more  clearly  the  various  elements 
which  comprise  the  cumulative  damage  model.  Figure  4.1  Is  a  block  diagram 
showing  the  infrastructure  of  the  methodology.  It  shows  that  there  are 
five  levels  of  development  in  the  model  construction.  Namely,  these  are 
(1)  numerical,  (2)  experimental,  (3)  theorization,  (4)  simulation,  and 
(5)  prediction/correlation. 

Numerical  Development.  The  major  objectives  In  the  numerical  development 
are  (1)  to  accurately  and  efficiently  compute  the  three-dimensional  stress 
field  in  a  given  laminate  with  or  without  a  prescribed  macrocrack  in  It, 
and  (2)  to  compute  the  crack  front  (or  crack-tip  in  one-dlmenslonal  crack 
propagation)  energy  release  rate  if  the  laminate  does  contain  a  prescribed 
macrocrack. 

Depending  on  the  nature  of  the  prescribed  problem,  a  two-dimensional 
and/or  three-dimensional  solution  method  must  be  used.  In  the  present 
research,  the  finite-element  method  Is  adapted  because  of  its  versatility 
In  handling  complicated  laminate  geometry,  as  well  as  loading  combinations. 

When  multiple  failures  are  to  be  simulated,  random  material  flaws  must 
be  represented  by  using  a  random  number  generator,  which  can  be  easily 
incorporated  in  the  crack  simulation  model. 

Experimental  Development.  Experiments  are  necessary  to  determine  the  per¬ 
tinent  material  properties  required  as  Input  data  for  the  analytical  model. 
Four  types  of  properties  are  required  which  must  be  obtained  from  indepen¬ 
dent  tests  using  the  basic  UD  material.  These  are  (1)  ply  deformation 
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moduli  (2)  fracture  toughness  for  matrix-dominant  crack  propa¬ 

gation  (G^)>  (3)  effective  flaw  distribution  parameters  for  multiple 
failures  In  laminates,  and  (4)  fatigue  growth  parameters  for  matrix- 
dominant  crack  growth  rate  under  fatigue  load.  All,  except  (1),  are 
still  In  their  early  stage  of  development.  But,  Increased  Information 
is  becoming  available. 

Theorization  Development.  For  static  crack  growth,  the  energy  release 
rate  criterion  Is  used;  for  fatigue  growth,  the  power  law  based  on  energy 
release  rate  for  the  crack  growth  rate  is  used;  and  for  growth  under  ran¬ 
dom  loading,  the  concept  of  constant  damage  lines  is  used.  Again,  all  are 
In  their  early  stage  of  development;  more  research  Is  needed  In  order  to 
settle  some  of  the  basic  questions  In  their  use. 

Simulation  Development.  This  is  the  actual  modeling  of  the  crack  growth 
process,  which  puts  together  the  use  of  the  computational  tool,  the  experi¬ 
mental  Input  data  for  material  properties,  and  the  basic  theory  for  crack 
growth.  The  model(s)  thus  produce  the  desired  output  Information  for 
which  the  reliability  of  the  laminate  Is  to  be  assessed. 

Prediction/Correlation.  The  predicted  outcome  Is  compared  by  experiment. 
The  comparison  serves  two  purposes  (1)  to  re-adjust  the  simulation  and/or 
theorization  developments  and  (2)  to  verify  the  applicability  of  the 
model (s)  and  establish  the  bounds  In  which  Is  Is  applicable. 
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APPENDIX  A 


COMPUTATION  OF  ENERGY  RELEASE  RATE  RETENTION  FACTOR  R(Sl 


Consider  the  Inset  In  Figure  A-1.  Assume  that 

there  exists  already  a  transverse  crack  which  Is  to  the  right  of  a  flaw  having  a 
size  of  2a.  The  effect  of  the  crack  on  the  energy  release  rate  behavior  of 
the  flaw  Is  a  reduction  of  the  energy  release  rate,  depending  on  the  distance 
S  between  the  two. 

Figure  A-1  shows  a  family  of  curves  for  the  energy  release  rate  coefficient 
C^  (as  a  function  of  flaw  size  a/t)  at  different  values  of  S.  It  Is  seen  that 
a  flaw  located  at  S  >  ISt,  nearly  lOOZ  of  the  energy  release  rate  Is  retained. 
This  Implies  that  the  existing  crack  no  longer  has  any  effect  on  the  flaw,  or 
that  the  flaw  Is  outside  the  shear-lag  zone  of  the  existing  crack. 

In  fact,  for  a  flaw  of  any  size,  the  influencing  factor  depends  only  on 
the  value  of  S.  Thus,  from  the  family  of  curves  such  as  shown  in  Figure  A-1, 
an  energy  release  rate  retention  factor  R(S)  may  be  determined.  For  the 
[02/902]^  and  Instance,  R(S)  curves  are  displayed  in 

Figure  A-2  as  plotted  against  S/t  (t  being  the  thickness  of  one  ply). 

The  R(S)  curves  for  [®2^^®n^s’  ®  fact  be  reduced  to  a 

single  curve  when  S  Is  normalized  by  nt,  see  Figure  A-3.  And,  this  single 
curve  can  be  fitted  by  an  exponential  function 


R(S/nt)  -  1  - 


1  <  S/nt  ^  8 
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Figure  A-3  Energy  Release  Rate  RetenUon  Factor  R(S)  As  A  fu  nctlon  of  Normalized  Distance  s/nt. 


APPENDIX  B.  TRANSVERSE  il^RACKING  -  STATIC  LOADING 


The  basic  material  used  In  the  experiment  was  AS-3501-06.  But, 
the  test  specimens  used  in  Phase-I  were  fabricated  (autoclave  curing) 
at  NADC  (Warminster,  PA);  while  the  test  specimens  used  in  Phase-II  were 
fabricated  at  the  University  of  Delaware.  The  basic  material  data 
between  the  two  batches  of  specimens  exhibited  a  slight  difference,  but 
enough  to  Indicate  that  they  must  be  treated  as  two  different  material 
systems  in  our  modelling  correlation. 

For  the  transverse  cracking  simulation,  two  sets  of  elastic  moduli 

were  used.  These  are  summarized  in  a  comparative  table  in  Appendix  H. 

For  mode-I  cracking,  the  critical  energy  release  rate  G^  was  set  at 

2 

G^  =  1.3  Ib/in  (228  J/m  )  for  both  batches  of  specimens.  The  effective 
flaw  size  distribution  f (a)  for  batch  1  and  batch  2  were  both  assumed  in 
the  form  of  normal  distribution.  But  their  individual  mean  and  standard 
deviation  were  slightly  d1f''  '-.nt; 


Batch  1 : 


Batch 


0.0032"  (0,08  mm) 
i).0008"  (0.02  mm) 

0.0026"  (0.007  mm) 
0.00042"  (0.011  mm) 


The  experimental  data  for  transverse  crack  growth  under  static  loading 

of  batch  1  specimens  are  shown  in  Figures  B-1  and  B-2  for,  respectively, 

[0  /90  ]  and  [0„/90.]  laminates.  The  Monte-Carlo  simulations  for  these 
2  2s  23s 

two  laminates  are  shown  in  Figures  B-3  and  B-4,  respectively. 
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The  experimental  data  for  transverse  crack  growth  under  static 
loading  of  batch  2  specimens  are  shown  in  Figures  B-5  and  B-6  for,  respec¬ 
tively,  corresponding  Monte-Carlo 

simulations  are  depicted  in  Figures  B-7  and  B-8. 
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Figure  B-1  Experimental  Transverse  Crack  Density  Versus  Applied  Laminate  Stress 
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Figure  B-3  Simulated  Transverse  Crack  Density  Versus  Applied  Laminate  Stress 
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Simulated  Transverse  Crack  Density  Versus  Applied  Laminate  Stress 
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Figure  3-5  Experimental  Transverse  Crack  Density  Versus  Applied  Laminate  Stress 
for  [02/902]g  Laminates. 
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Figure  B-8  Simulated  Transverse  Crack  Density  Versus  Applied  Laminate  Stress 
for  [02/902]s  Laminates. 


APPENDIX  C.  TRANSVERSE  CRACKING  -  FATIGUE  LOADING 


In  Phase-I,  transverse  cracking  under  constant  amplitude  fatigue 

loading  was  studied  for  [0„/90.]  and  [0„/90.]  laminates.  Experiments 

Z  Z  s  Z  ^  S 

were  performed  under  amplitudes  =  26  ksi  (179  mpa) ,  30  ksi  (207  mpa) , 
and  38  ksi  (262  mpa)  for  under  amplitudes  0^  =  38  ksi 

(262  mpa),  A3  ksi  (297  mpa),  and  55  ksi  (365  mpa)  for 
experimental  data  are  shown  graphically  in  Figure  C-1  for  t®2^^^2^s 
Figure  C-2  for  t^2^^^3^s‘ 

Simulation  results  using  the  batch  1  material  constants  and  the  crack 
growth  equation 


da 

dN 


G(0  ,  a)  - 

0l[ - f 


(C-1) 


with 

a  =  O.OOAt,  p  =  20  (C-2) 

are  shown,  respectively,  in  Figures  C-3  and  C-A. 

In  Phase-II,  the  same  laminates  were  tested  as  in  Phase-I,  except 
that  the  batch  2  specimens  were  characterized  with  a  set  of  different 
material  properties  as  discussed  earlier. 


Experiments  under  constant  amplitude  fatigue  loading  were  performed 

under  =  35  ksi  (2A1  mpa),  A5  ksi  (310  mpa),  60  ksi  (A13  mpa),  68  ksi 

(A68  mpa),  and  76  ksi  (52A  mpa)  for  under  25  ksi  (172  mpa), 

35  ksi  (2A1  mpa),  A5  ksi  (310  mpa),  50  ksi  (3A5  mpa),  and  55  ksi  (379  mpa) 

for  [0-/90  1  . 

2  3  s 

Figures  C-5  through  C-8  show  graphically  the  experimental  data  ob¬ 
tained  for  [0-/90_]  .  ;.'nte  that  data  for  CT  =  35  ksi  was  not  shown.  The 

2  z  s  t 
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reason  Is  that  no  transverse  cracks  were  found  under  this  loading  for 


up  to  one  million  cycles. 

For  the  same  reason,  data  for  =  25  ksi  and  =  35  ksi  in  the 

[0„/90_]  tests  were  also  not  shown.  Data  obtained  under  the  other  three 
2  is 

amplitudes  are  shown  in  Figures  C-9  through  C-11. 

Fatigue  simulations  of  the  cases  shown  in  Figures  C-5  through  C-11 
were  conducted  using  the  material  constants  obtained  for  the  batch  2 
specimens.  The  fatigue  crack  growth  equation  remained  as  in  Equation  C-1; 
and  the  constants  Ot,  p  also  remained  as  in  Equation  C-2. 


The  corresponding  simulated  results  are  shown  in  Figures  C-12  through 


C-18. 
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Figure  c-1  Simulated  Transverse  Crack  Density  Versus  Fatigue  Cycle  N.  £o  /90  ] 
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Figure  C-2  Experimental  Transverse  Crack  Density  Versus  Fatigue  Cycl 
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Figure  C-3  Experimental  Transverse  Crack  Density  Versus  Fatigue  Cycle  N.  [o^/90^J 
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.EXPERinENT 


Experimental  Transverse  Crack  Density  Versus  Fatigue  Cycles 
for  [0„/90-]  Laminates.  Maximum  Stress  of  O  =45  KSI. 


EXPERIMENT 


Figure  C-6  Experimental  Transverse  Crack  Density  Versus  Fatigue  Cycles 
for  [0-/90„]  Laminates.  Maximum  Stress  of  a  *  60  KSI. 
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Figure  C-7  Experimental  Transverse  Crack  Density  Versus  Fatigue  Cycles 
for  [0_/90-]  Laminates'.  Maximum  Stress  of  a  -68  KSI. 
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Figure  C-8  Experimental  Transverse  Crack  Density  Versus  Fatigue  Cycles 
for  [0-/90_]  Laminates.  Maximum  Stress  of  a  =76  KSI. 
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Figure  C-9  Experimental  Transverse  Crack  Density  Versus  Fatigue  Cycles 
for  fO  /90  ]  Laminates.  Maximum  Stress  of  a  =  45  KSI. 


Figure  C-10  Experimental  Transverse  Crack  Density  Versus  Fatigue  Cycle 
for  [0_/90  ]  Laminates.  Maximum  Stress  of  Cf  =50  KSI. 
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igure  C-11  Experimental  Transverse  Crack  Density  Versus  Fatigue  Cycles 
for  [0-/90-]  Laminates.  Maximum  Stress  of  a  =55  KSl. 
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Figure  C-12  Simulated  Transverse  Crack  Density  Versus  Fatigue  Cycles 
for  [0^/90- ]  Laminates.  Maximum  Stress  of  o  =45  KSI. 
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Figure  C-13  Simulated  Transverse  Crack  Density  versus  Fatigue  Cycles 
for  [0„/90„]  Laminates.  Maximum  Stress  of  a  =60  KSI. 


Figure  C-14  Simulated  Transverse  Crack  Density  Versus  Fatigue  Cycles 
for  [Oo/^O-]  Laminates.  Maximum  Stress  of  o  “68  KSI. 
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Figure  C-15  Simulated  Transverse  Crack  Density  Versus  Fatigue  Cycles 
for  ^®“inates.  Maximum  Stress  of  a  =76  KSI. 
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Figure  C-16  Simulated  Transverse  Crack  Density  Versus  Fatigue  Cycles 
or  Maximum  Stress  of  a  =45  KSI. 


CUG/4»]  Airs 

(S) 

r-H 


CNI/#]  AiESNBQ  >IGyd3 


-  129 


Figure  C-18  Simulated  Transverse  Crack  Density  Versus  Fatigue 
for  [0„/90_]  Laminates.  Maximum  Stress  of  0  “55  KSI. 


APPENDIX  D.  TRANSVERSE  CRACKING  -  CUMULATIVE  DAMAGE 

Dl,  Phase-I  Data 

The  simulation  results  of  transverse  crack  growth  under  static  and 
constant  amplitude  fatigue  loading  were  presented  in  Appendices  B  and  C. 
There,  the  material  effective  flaw  size  distribution  in  the  90“-layer  was 
assumed  in  the  form  of  a  normal  distribution,  with  the  parameters  (mean 
and  standard  deviation) : 

U  =  0.0032"  (0.08  mm) 

=  0.0008"  (0.02  mm) 

For  fatigue  growth,  the  parameters  o  and  p  in  the  growth  rate 
equation  were  chosen  as 

a  =  0.004t,  P  “  20. 

Using  these  parameter  values  in  the  simulation  routines  for  fatigue 

growth,  a  set  of  "constant  damage"  lines  was  generated  in  the  o  -  N  space 

for  each  of  the  [0„/90„]  and  the  [0_/90_]  laminates.  These  two  sets  of 
/  /  s  z  J  s 

"constant  damage"  lines  are  shown  in  Figures  D-1  and  D-2. 

For  cumulative  damage  growth  under  variable  amplitude  fatigue,  both 

the  [0-/90- ]  and  [0-/90_]  laminates  were  tested  in  the  experiment.  The 
2  2s  23s 

following  test  cases  were  carried  out: 

[0-/90- 1  :  0-  =  60  ksi  (414  mpa)  for  10^  cycles 

2  2s  f 

=  80  ksi  (552  mpa)  iol  10^  cycles 

2 

Oj  =  80  ksi  (552  mpa)  for  10  cycles 

60  ksi  (414  mpa)  for  4  x  10^  cycles 


[0„/90_]  ; 
2  2  s 


[0-/90  ]  : 
2  3  s 


50 

ksi 

(345 

mpa) 

for 

10^ 

cycles 

70 

ksi 

(483 

mpa) 

for 

10^ 

cycles 

50 

ksi 

(345 

mpa) 

for 

10^ 

cycles 

26 

ksi 

(179 

mpa) 

for 

lo'^ 

cycles 

30 

ksi 

(207 

mpa) 

for 

3  X 

4 

10  cycles 

26 

ksi 

(179 

mpa) 

for 

10^ 

cycles 

15 

ksi 

(104 

mpa) 

for 

10^ 

cycles 

50 

ksi 

(345 

mpa) 

for 

10^ 

cycles 

80 

ksi 

(552 

mpa) 

for 

10  cycles 

40 

ksi 

(276 

mpa) 

for 

10^ 

cycles 

Transverse  cracks  per  unit  length  specimen  were  determined  after  each 
load-step.  These  are  then  compared  to  the  same,  which  are  predicted  by 
the  cumulative  damage  model  using  the  constant  damage  lines  shown  in 
Figures  D-1  and  D-2.  The  comparison  is  illustrated  in  Tables  D-1  and  D-2. 
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TABLE  D-1.  CUMULATIVE  DAMAGE  DATA  FOR  PHASE-I 


r 

( 

[ 

■ 


Laminate 


Damage ,  Ck/ In 


Fatigue  Loads 

Cycles  at  Load 

Experiment 

Prediction 

60  ksi 

10^ 

33 

27 

80  ksi 

10^ 

36 

37 

60  ksi 

10^ 

34 

27 

80  ksi 

10^ 

36 

37 

60  ksi 

10^ 

34 

27 

80  ksi 

10^ 

37 

37 

80  ksi 

10^ 

29 

29 

60  ksi 

10^ 

36 

31 

80  ksi 

10^ 

29 

29 

60  ksi 

10^ 

34 

31 

80  ksi 

10^ 

31 

29 

60  ksi 

10^ 

37 

31 

50  ksi 

10^ 

19 

13 

70  ksi 

10^ 

31 

28 

50  ksi 

10^ 

32 

29 

50  ksi 

10^ 

11 

13 

70  ksi 

10^ 

28 

28 

50  ksi 

10^ 

30 

29 

50  ksi 

10^ 

16 

13 

70  ksi 

10^ 

28 

28 

50  ksi 

10^ 

30 

29 

► 


I. 


I 
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Laminate 


[02/903]^ 


TABLE  D-2.  CUMULATIVE  DAMAGE  DATA  FOR  PHASE-I 


Fatigue  Loads 


Cycles  at  Load 


4 


26 

ksl 

10 

30 

ksi 

3 

X  10 

26 

ksl 

10^ 

26 

ksi 

10' 

30 

ksi 

3 

X  10' 

26 

ksl 

10^ 

26 

ksi 

lo"* 

30 

ksi 

3 

X  10' 

26 

ksi 

10^ 

26 

ksi 

10^ 

30 

ksl 

3 

X  10' 

26 

ksi 

10^ 

15 

ksl 

10^ 

50 

ksl 

10^ 

15 

ksl 

10^ 

50 

ksl 

10^ 

80 

ksl 

10 

40 

ksl 

10^ 

80 

ksl 

10 

40 

ksl 

10^ 

Damage.  Ck/ln _ 

Experiment  Prediction 


0 


D2.  Phase-II  Data. 


» 


The  material  effective  flaw  size  distribution  in  the  90°-layer  for 
the  Phase-II  study  was  also  assumed  in  the  form  of  a  normal  distribution. 
But,  the  parameters  were  chosen  as  follows: 

n  =  0.0026"  (0.07  mm) 

'^a 

a  =  0.00042"  (0.011  mm) 

a 

All  other  values  used  in  the  simulation  were  the  same  as  in  Phase-I. 

The  constant  damage  lines  for  [0„/90„]  and  [0„/90,]  laminates  wc.re  also 

z  z  s  Z  J  s 

generated  by  the  same  simulation  routine.  These  are  shown  in  Figures  D-3 
and  D-4,  respectively. 

Cumulative  damage  tests  on  these  two  laminates  were  carried  out  for 
the  following  cases: 

[02/902]^:  =  60  ksi  (345  mpa)  for  10^  cycles 

80  ksi  (552  mpa)  for  10^  cycles 

2 

=  80  ksi  (552  mpa)  for  10  cycles 

2 

60  ksi  (414  mpa)  for  10  cycles 

3 

0^  =  50  ksi  (345  mpa)  for  10  cycles 

3 

70  ksi  (483  mpa)  for  10  cycles 

50  ksi  (345  mpa)  for  10^  cycles 

L 

[02/90^]^:  cjj  “  26  ksi  (179  mpa)  for  10  cycles 

4 

30  ksi  (207  mpa)  for  2  x  10  cycles 

26  ksi  (179  mpa)  for  10^  cycles 

=  15  ksi  (104  mpa)  for  10^  cycles 

50  ksi  (345  mpa)  for  10^  cycles 
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AO  ksi  (276  mpa)  for  10^  cycles 

Cumulative  damage  In  tenns  of  transverse  crock  growth  (cracks  per 
inch)  are  determined  both  from  experiment  and  from  the  cumulative  damage 
model  (using  the  constant  damage  lines). 

The  comparison  between  them  is  tabulated  in  Table  D-3  for 
and  in  Table  D-A  for  [0„/90.]  . 


TABLE  D-3.  CUMULATIVE  DAMAGE  DATA  FOR  PHASE-II 


Damage.  Ck/in 


Laminate 

Fatigue  Loads 

Cycles  at  Load 

Experiment 

Prediction 

60  ksi 

10^ 

22 

21 

80  ksi 

10^ 

33 

33 

60  ksi 

10^ 

13 

21 

80  ksi 

10^ 

34 

30 

(Failed  before 

10^  cycles) 

80  ksi  10^ 

60  ksi  10^ 


80  ksi  10^ 

60  ksi  10^ 


50  ksi 
70  ksi 
50  ksi 


10 

•3 

10" 

10' 


10' 

10^ 

10^ 


50  ksi 
70  ksi 
50  ksi 
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CYCLES  L0GCN3 

Figure  D-2  Constant  Damage  Lines  for  Batch  1  [0,/90-]  Laminates. 
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Figure  D-4  Constant  Damage  Lines  for  Batch  2  [0  /90  ]  Laminates. 


APPENDIX  E.  FREE  EDGE  DELAMINATION  -  STATIC  LOADING 


El .  Phase-I  Data. 

Free  edge  delamination  data  obtained  under  static  loading  in  Phase-I 
study  has  already  been  presented  in  Section  2.6.  There,  the  laminates 
[+25/90^]^,  n  =  1,  2,  3,  were  tested  and  simulated.  It  was  noted  that 
the  cases  of  n  =  2  and  3  involved  transverse  crack/delamination  interactions. 
Consequently,  a  vigorous  delamination  simulation  for  these  two  cases  was  not 
available.  However,  the  onset  of  delamination  loads  for  the  latter,  were 
determined  approximately  (see  discussions  in  Section  2.6  and  Section  2.9). 

In  order  to  avoid  duplication,  these  results  will  not  be  Included  here. 

E2 .  Phase-II  Data. 

In  Phase-II,  free  edge  delamination  study  was  conducted  also  using 
the  laminates  [+25/90^]^,  n  =  1,  2,  3.  For  the  cases  of  n  =  2  and  3, 
similar  transverse  crack/delamiaation  interaction  was  also  encountered.  In 
the  simulation,  a  new  set  of  energy  release  rate  curves  was  generated, 
using  the  new  material  moduli  for  the  batch  2  material.  These  curves  are 
shown  in  Figures  E-1  and  E-2. 

By  means  of  the  energy  release  rate  curves  and  the  Griffith  crack 

growth  criterion,  the  predicted  onset  loads  for  free  edge  delamination  are 

2 

compared  to  experiment  as  listed  in  the  following,  (G^  =  1.0  in-lb/in  ): 


Onset  Laminate 

Strain,  e 

X 

Laminate 

Experiment 

Prediction 

[+25/90] 

—  s 

5.78  X  10"^ 

5.68  X  10"^ 

[+25/90J 
—  2  s 

5.52  X  10"^ 

5.68  X  10“^ 

[+25/90J 
“  J  s 

5.07  X  10”^ 

5.22  X  10"^ 

For  details  on  the  calculated  results,  refer  to  Section  2.6  for 
[+25/90]^.  and  Section  2.9  for  [+25/90,]  and  [+25/90  ] 

s  —  3 ■'g 

Some  graphical  data  pert.lnlag  to  the  above  may  ba  found  l„  Appendix  H 


20X10"® 


Curves  For  C±25/90^]s  Laminates 
Batch  2  ,Ct 


Figure  E-2  Energy  Release  Rate  Coefficients  for  Delamination 
Under  Thermal  Load. 


APPENDIX  F.  FREE  EDGE  DELAMINATION  -  FATIGUE  LOADING 

Fl.  Phase-I  Data. 

Again,  the  delamination  growth  data  under  constant  amplitude  fatigue 

loading  obtained  in  Phase-I  has  already  been  presented  in  Section  2,7. 

Recall  that  only  the  laminate  [+25/90]^  was  tested  and  simulated.  The 

simulated  results  were  based  on  the  crack  growth  in  the  25/90  Interface  and 
by  means  of  the  equation 

§  -  5[G(o^.  a)/GjP  (F-1) 

where  the  parameters  were  selected  as 

a  =  0.04  inch;  p  =  8.3.  (F-2) 

2 

Given  the  G(a^,  a)  curve,  and  setting  G^  =  1.3  in-lb/in  ,  the  crack 
growth  equation  can  be  integrated  numerically  to  obtain  the  quantity  a  (size 
of  delamination)  in  terms  of  0^  and  N.  All  relevant  results  are  found  in 
Section  2.7. 

F2 .  Phase-II  Data. 

The  [+25/90]  laminates  made  of  batch  2  material  were  tested  and  simu- 
—  s 

lated  for  free  edge  delamination  under  the  following  fatigue  loads: 

Oj  *  29  ksl  (200  mpa) ;  33  ksi  (278  mpa) ; 

40  ksl  (276  mpa);  45  ksi  (310  mpa); 

50  ksi  (345  mpa);  55  ksi  (379  mpa). 

The  tested  and  simulated  results  for  the  above  are  shown  in  Figures 
F-1  to  F-6.  All  simulations  were  based  on  the  growth  equation  (F-1),  with 
the  same  parameter  values  in  (F-2) . 
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Percent  Delomination 


Figure  F-1  Percent  Delamlnatlon^vs.  Fatigue  Cycles  for 
(+25/90 Laminates,  o  •  29  KSI. 
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Percent  Delomination 


Cycles 


figure  t“"2  Percent  DelamLnation  vy.  Fatigue  Cycles  for 
[+25/90]^  Laminates,  o  =  33  KSI. 
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Percent  Delominotion 


Cycles 


Figure  F-3  Percent  Delamination  vs.  Fatigue  Cycles  for  Laminates 
a  =  40  KSI. 


Percent  Delominotion 


Cycles 


Figure  F-A  Percent  Delamination  vs.  Fatigue  Cycles  for 
[+25/90]^  Laminates,  o  =  A5  KSI. 
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C±25/903s  Lominote  Delominotion 
Under  Fatigue  Loading,  Batch  2 


I  •  Data 
r  ''Theory 


•  •• 


p 


Percent  Delamination 


C±25/90]s  Laminate  Delamination 
Under  Fatigue  Loading,  Botch  2 


9 


Figure  F-6  Percent  Delamination  vs.  Fatigue  Cycles  ior 
[+25/90]^  Laminates,  a  =  55  KSI. 
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APPENDIX  G.  FREE  EDGE  DELAMINATION  -  CUMULATIVE  DAMAGE 


Gl.  Phase-I  Data. 

Some  of  the  free  edge  delamlnatlon  data  under  variable  amplitude 
fatigue  loading  obtained  In  Phase-I  was  used  In  the  discussions  In 
Section  2.8,  For  completeness,  all  test  results  and  their  simulated 
counterparts  are  graphically  displayed  In  Figures  G-1  to  G-5. 

Note  that  In  each  of  the  figures,  the  test  results  and  the  simulated 
results  are  displayed  In  the  a  -  N  space,  which  is  spanned  by  a  set  of 
"constant  damage  lines."  The  constant  damage  lines  were  generated  by 
the  crack  growth  equation  (F-1)  along  with  the  parameters  in  (F-2) .  Test 
results  (delamlnatlon  in  %  of  specimen  half-width)  in  each  of  these 
figures  are  Indicated  by  the  number  (in  %)  in  parentheses. 
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Cumulative  Damage  Curves  For  C±25/903 
Laminates,  Batch  1 


Figure  G-2  Cumulative  Damage  Tests  for  [+25/90 Laminates. 


Cumulative  Damage  Curves  For  C±25/90] 
Lamlnofes,  Batch  I 


Cumulative  Damage.  Test  for  [+25/90 3^  Laminates. 


G2.  Phase-II  Data. 


Similar  comparative  results  obtained  in  Phase-II  are  shown  in 
Figures  G-6  through  G-10.  Here,  again,  are  five  study  cases  for  the 
cumulative  damage  under  variable  amplitude  fatigue  loadings. 


Cumulative  Damage  Curves  For  C±25/903 
Laminates,  Batch  2 

20  %  \  90% 


Cumulative  Damage  Tests  for  1+2J/90]  Laminates 


Cumulative  Damage  Curves  For  C±25/903 

Laminates,  Batch  2 
^  * 


Cumulative  Damage  Tests  for  [+25/90 Laminates 


APPENDIX  H.  BASIC  MATERIAL  CHARACTERIZATION 


Hi.  Material  Property  Characterization  in  Phase-II. 

A  new  batch  of  material  (AS-3501-06)  was  used  in  the  experiment  during 
Phase-Il.  Material  properties  of  the  new  batch  were  determined  and  compared 
with  those  obtained  during  Phase-1. 

Some  of  the  important  material  properties  for  the  U.D.  and  other 
laminates  are  shown  in  Tables  H-1,  H-2,  and  H-3.  Generally  a  +10%  differ¬ 
ence  is  noted.  But,  this  difference  is  considered  to  be  within  the  experi¬ 
mental  scatter.  So,  basically,  the  new  batch  of  material  is  considered 
similar  to  the  old  batch  of  material. 

Figure  H-1  illustrates  the  tensile  strength  property  scatter  of  [Og] 
laminate  between  the  old  and  new  batches  of  materials.  In  this  case,  the 
new  batch  of  material  exhibited  a  great  scatter  in  the  tensile  strength 
property. 

A  comparison  of  the  transverse  crack  growth  curves  (number  of  cracks 
versus  applied  tensile  load)  of  t02'^^®2^s  batches 

is  shown  in  Figure  H-2;  a  similar  comparison  for  t02/90g]^  is  shown  in 
Figure  H-3.  The  present  batch  showed  a  stronger  resistance  to  transverse 
cracking.  This  is  due  possibly  to  the  slightly  stiff er  E  property  of  the 

ij 

0®-ply,  see  Table  H-1. 

The  free  edge  Induced  delamlnatlon  growth  curves  in  [+25/90  ]  ,  n  =  1, 

—  ns 

2,  3  are  shown,  respectively,  in  Figures  H-4,  H-5,  and  H-6.  Again,  the 
difference  between  the  new  and  the  old  batches  can  be  noted.  Generally, 
the  new  batch  of  specimens  developed  edge  delamlnatlon  at  a  lower  load 
and  the  growth  of  it  also  showed  a  faster  rate.  The  reason  for  the  differ¬ 
ence  is  not  immediately  clear,  however. 
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■eH-1  Compaxlson  of  the  Tensile  Strength  Distributions  of  the  New  and  the  Old  Batches 
of  Materials. 


batch 


A  Compaxison  of  the  Transverse  Crack  Growth  Curves  Between  the  Hew  and  the  Old 
Batches  of  Materials.  r0o/90_1 


Previous  batch 


Figure  H-6  A  Comparison  of  the  Free  Edge  Delamination  Growth  Curves  Between  the  New 
and  the  Old  Batches  of  Materials,  C+25/90^3_ 


